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CHAPTER  1 
INTRODUCTION 


The  Naval  Surface  Warfare  Center,  Dahlgren  Division  (NSWCDD)  has  been  one  of  the 
prime  forces  behind  the  resonant  scattering  theory  (RST)  venture  during  the  last  two  decades. 
NSWCDD  patronage  and  sponsorship  led  to  numerous  publications  in  archive  journals’"^  and 
In  book  chapters.^^"^  At  least  three  U.  S.  Patents^®^  were  granted  as  a  result  of  research  and 
development  (R&D)  efforts  on  the  RST  and  related  topics.  Notable  contributions  are  attributed 
to  G.  Gaunaurd  (R42)  and  collaborators  D.  Brill  (The  United  States  Naval  Academy)  and  H. 
Uberall  (Catholic  University  of  America).  Experimental  verification  and  development  by  the 
authors  of  this  report  started  several  years  ago  in  the  Underwater  Systems  Department  (U).  The 
Instrumentation  and  Electronics  Branch  (U42)  portion  of  that  department  where  this  work  was 
done  was  reorganized  to  the  Panama  City,  Coastal  System  Station,  Code  2620,  Mine  Systems 
Branch. 

The  concept  behind  active  sonar  and  the  RST  is  deceptively  simple  and  has  been  known 
since  the  advent  of  classical  physics.  When  a  load  is  applied  to  a  structure,  the  structure 
reacts  accordingly.  According  to  Newton’s  laws,  objects  at  rest  tend  to  remain  at  rest  and  if  in 
motion,  tend  to  stay  in  motion.  Relating  this  to  a  structure,  when  a  sinusoidal  load  Is  applied  to 
a  structure,  most  of  the  energy  will  be  reflected  until  steady  motion  is  achieved.  The  body  will 
then  vibrate  at  all  the  resonant  modes  that  fall  within  the  power  spectrum  of  the  load.  The 
energy  associated  with  these  modes  will  be  absorbed  and  stored  in  the  structures  and  the 
remainder  of  the  energy  will  be  reflected. 

The  main  thrust  of  this  research  was  to  experimentally  determine  the  resonances  in 
cylindrical  like  structures.  When  such  a  structure  Is  insonified  with  a  single  frequency 
sinusoidal  acoustic  energy  there  will  be  some  deformation  of  the  structure  and  internal  waves 
will  propagate  through  it.  Some  waves  will  travel  from  one  edge,  through  the  center,  to  the 
other  edge,  and  be  reflected  back.  Others  will  travel  around  the  structure  circumferentially.  If 
the  driving  frequency,  the  radial  waves,  and  the  circumferential  waves  are  in  phase  at  the 
driving  point,  resonance  will  occur.  The  circumferential  waves  are  usually  described  by 
sinusoidal  functions  while  the  radial  waves  are  described  by  Bessel  and  Hankel  functions.  If 
the  Incident  frequency  increases,  higher  resonance  frequencies  will  occur  in  both  waves,  thus 
Indices  must  be  applied.  For  the  circumferential  waves  the  index  of  n  will  be  used,  and  for  the 
radial  waves  /.  Since  the  necessary  condition  for  resonance  is  that  both  the  circumferential  and 
the  radial  be  excited,  resonances  are  classified  by  their  pairs  of  indices,  (n,/).  When  the  n  index 
is  determined  experimentally,  the  number  of  standing  waves  developed  around  the 


1-1 


css  TR  456-93 


circumference  of  the  cylinder  is  known.  Basically  the  resonance  of  the  structure  is  dependent 
on  the  size,  shape  and  material  it  is  manufactured  from.  A  review  of  theory  reveals  that  the 
velocity  of  sound  in  material  is  the  square  root  of  the  ratio  of  one  of  the  material  moduli  and  the 
materi^  density  c  =  Vl/p.  The  wave  number  of  the  material  is  (k  =  27tf/c)  and  for  a  cylinder,  the 
radius  (a)  is  the  important  size,  each  resonant  frequency  will  provide  a  data  point  in  the  ka 
domain  and  a  set  of  these  resonances  will  identify  the  target.  Knowing  that  targets  are  only 
made  out  of  a  limited  class  of  materials  and  a  limited  number  of  sizes,  a  library  of  resonances 
for  possible  targets  can  be  determined  theoretically  and  compared  to  the  measured  values. 

Stated  another  way  the  dynamic  response  at  any  point  of  the  structure  is  the  sum  of  the 
individual  modal  responses.  These  responses,  fashioned  after  the  modal  power  partitions  from 
the  total  input,  are  governed  by  the  so-called  frequency  equation  which  is  a  function  of  the 
structural  parameters  and  the  input  frequency.  In  the  cases  where  one  input  frequency 
matches  one  of  the  structural  natural  frequencies,  the  mode  associated  with  this  frequency  is 
vigorously  excited  to  result  in  a  resonance  vibration.  At  resonance,  the  majority  of  the  input 
power  is  directed  to  this  excited  mode  to  sustain  its  vigorous  motion  at  the  expense  of  the 
others.  Should  the  input  last  for  a  while,  the  structure  will  enter  into  a  steady  motion.  Energy 
continues  to  enter  the  structure  only  to  be  reradiated  due  to  the  motion  of  the  body,  be  stored 
as  kinetic  energy  or  be  converted  into  heat  through  a  frictional  mechanism  internally  and  is 
released  into  the  surroundings.  At  the  cessation  of  the  input,  the  structure  will  revert  to  free 
vibration  and  the  stored  energy  will  eventually  be  dissipated. 

If  the  assumption  is  made  that  the  structure  in  water  is  lightly  damped,  this  will  be 
discussed  later,  the  behavior  of  structures  underwater  is  the  same  as  that  described  in  the  last 
paragraphs.  When  a  wave  insonifies  a  target,  energy  of  the  wave  input  penetrates  into  the 
target  and  sets  its  internal  elements  in  motion.  The  balance  of  the  insonification  is  reflected. 

The  reflection  yields  the  bistatic  crosssection  which  can  be  measured  by  a  hydrophone  slowly 
revolving  around  the  structure  at  the  far  field.  Normally,  the  bistatic  crosssection  has  a  very 
complicated  pattern  because  it  is  a  result  of  a  conglomeration  of  many  partial  waves  scattered 
by  the  target  at  broadside.  It  also  has  a  large  magnitude  since  most  energy  is  reflected  and 
little  penetrates  into  the  target.  At  a  resonant  state,  the  monochromatic  bistatic  crosssection  will 
display  a  much  simpler,  symmetrical  pattern  known  as  the  "rosetta."  This  rosetta  pattern  has  a 
simple  mathematical  description  so  an  amplitude  measured  at  the  maximum  value  suffices  to 
specify  the  amount  of  the  reradiated  power  and,  by  deduction,  the  amount  which  is  penetrated 
into  the  target.  The  magnitude  of  the  broadside  steady  state  reflection  is  expected  to  be  much 
smaller  at  resonance  since  most  of  the  input  power  is  penetrating  the  target. 

The  comparison  of  energy  returned  from  the  target  at  resonance  frequency  (when  it  is 
compliant)  and  off  resonance  (when  it  is  rigid)  gives  rise  to  the  idea  of  experimentally  locating 
the  unique  set  of  resonances  which  is  the  signature  of  the  structure.  Determining  the  location  of 
several  of  these  frequencies  is  sufficient  to  ascertain  the  signature  for  the  target  identification 
and  classification  purposes. 
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The  goal  of  the  investigation  reported  here  was  to  develop  a  methodology  for  the 
extraction  of  energy  data  from  the  return  at  resonance  and  off  resonance  states  such  that  the 
signature  of  a  target  underwater  can  be  measured.  This  measured  signature  contains  the 
resonance  frequencies  of  the  target  which  identifies  its  shape,  size,  and  composition  which  can 
be  used  for  its  identification  and  classification.  It  was  decided  that  the  simplest  targets  would 
be  used  because  experimental  results  can  be  verified  from  the  readily  available  predictions  of 
the  theory.  A  successful  completion  of  this  endeavor  will  validate  the  assumptions  and 
approaches  taken  in  the  RST  development  and  will  provide  a  solid  foundation  for  further 
research  directed  toward  securing  signatures  for  complicated  and  realistic  targets.  This  has 
important  and  practical  implications  since  the  more  complex  or  real  life  targets  are  impracticai 
to  model  analytically,  the  experimental  methodology  will  be  the  only  means  of  securing  the 
signatures. 

The  following  chapters  wili  provide  a  brief  description  of  the  theoretical  approach  to  the 
RST  solution  by  way  of  an  example,  followed  by  the  experimental  arrangements  and 
procedures  and  by  selected  experimental  results.  Recommendation  for  future  investigation  is 
also  presented. 
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CHAPTER  2 

RST  AND  ITS  FORMULATION 


The  equations  leading  to  the  RST  formulation  are  derived  from  standard  formulas  in  linear 
theories  of  elasticity  on  isotropic  solids  and  in  linear  acoustics  in  an  inviscid  fluid  media.  The 
two  sets  of  formulas  are  linked  by  boundary  conditions  at  the  interface,  namely,  the  pressure  at 
the  fluid  side  must  be  equal  to  the  normal  stress  in  the  solid  side  and  the  normal  displacement 
at  both  sides  must  be  equal.  In  addition,  since  an  inviscid  fluid  can  not  sustain  shear,  the  shear 
stress  at  the  interface  must  vanish.  Solving  this  boundary-value  problem  for  the  pressure 
distribution  in  the  fluid  medium  has  been  performed  for  many  geometrically  simple  targets  such 
as  infinitely  long,  elastic  cylinders,’®""**  spheres, “  cylindrical 
shells,’’^’’®‘^°'^’'^^  '^®’  “  hollow  spheres’^-*-®^  and  spheroids,®*  where  the  insonification  by  a  plane 
wave  is  at  a  direction  perpendicular  to  a  geometrical  axis  of  the  target.  For  an  illustration  of  the 
analytical  approach  taken  in  the  past,  the  solution  for  a  solid  cylinder,  infinite  in  length,  is 
presented  here  in  detail.  With  respect  to  Figure  2-1 ,  a  plane  wave  impinges  on  a  cylindrical 
target  of  radius,  a.  In  a  direction  perpendicular  to  its  axis  which  extends  to  Infinity  in  both 
directions. 
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FIGURE  2-1 .  INSONIFICATION  OF  A  CYLINDER 
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The  sound  pressure  of  interest,  p,  is  two-dimensional,  being  a  function  of  the  distance 
from  the  axis,  r,  and  the  angular  orientation  with  respect  to  the  direction  of  propagation  of  the 
plane  wave,  0.  By  this  definition,  the  backscattered  wave  (echo  or  reflected  wave)  is  at  an 
orientation  of  0  =  180“  or  it  radians.  In  the  fluid  side,  the  sound  pressure,  p(r,0,t).  is  governed  by 
the  wave  equation  in  cylindrical  coordinates: 


ar*  rdr  r*a0^ 


(2-1) 


where  c  is  the  speed  of  wave  propagation.  Assuming  a  harmonic  solution,  the  time  dependent 
term  in  p(r,0)e^"‘“*\  where  w  is  the  angular  velocity  of  the  incident  wave  vanishes  resulting  in  the 
Helmholtz  equation: 


£pl  i££i 


+kfp,  =  0 


(2-2) 


where  k,  =  <0/0,  is  the  wave  number  and  the  subscript  1  is  added  to  designate  the  fluid  medium. 
The  solution  of  Equation  (2-2)  is: 


P,(r,  0)  =  Po  i  (0"8„[J«(k|  r) + K Hi'>(k,  r)]  cosn 0  (2  -  3) 

n  =  0 


where  n  is  a  positive  integer  due  to  circumferential  periodicity;  e„  =  2  for  n  0  and  eo  =  1 .  In 
Equation  (2-3),  the  first  term  is  the  incident  plane  wave  of  amplitude  po  and  the  second  term, 

involving  the  (complex)  Hankel  function  of  the  first  kind,  gives  the  scattered  wave.  This  term 
consists  of  an  infinite  set  of  partial  waves,  each  of  amplitude  b„  whose  magnitude  can  be 

determined  from  given  boundary  conditions.  Once  b„  is  known,  the  sound  pressure  at  any  point 

in  medium  1  and  at  any  time  is  given  by:  p,(r,0)e^'‘“'\  As  an  example,  if  the  boundary  at  r  =  a  is 

rigid  or  immovable,  allowing  no  energy  penetration  into  the  solid  interior,  the  sole  condition 
which  remains  in  effect  is  the  vanishing  of  the  normal  displacement  at  the  boundary. 
Mathematically,  this  is  given  as  u,(a,0)  =  0  where  u,  is  the  radial  displacement.  Then  from 

Equation  (2-3),  since: 


u,(r,0) 


1  gPi 
p,o)^  5r 


(2-4) 


where  p,  is  the  density  of  the  fluid, 
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and.  at  the  boundary, 

'■,(«.e)  =  0=;^Xi‘4ji(k,a)+b.H<"'(k,a)]cosiie 

|./|  V]  CU  n  sO 


which  gives: 


lW  _  ■^nO^l  ft) 

"  Hi%.a) 


(2-5) 


The  primes  designate  derivatives  with  respect  to  the  argument  of  the  function.  Thus: 

00  J* 

piV.6)  =  Po  2(0X  *-  Hi'^Ocir)  cosnB  (2-6) 

«=o  Y  Hi,'>'(k,a)  J 

A  superscript,  (r) ,  is  added  to  designate  that  the  equation  is  pertaining  to  a  rigid  target.  If  the 
boundary  is  nonrigid,  hence  movable  as  in  the  case  for  what  is  known  as  an  elastic  target,  the 
movement  of  the  boundary  creates  a  stress  field  within  the  target.  The  governing  equation,  in 
terms  of  the  particle  displacement  vector,  u,  is: 

V"U+^V(V-5)— ^||=0  (2-7) 

where  X  and  ^  are  the  Lame's  constants  of  the  target  material  whose  density  is  p  and  c,  is  the 
speed  of  the  shear  wave,  (c, = ).  In  solid,  there  exists  additionally  a  dilatational  wave 
whose  speed  ca  is  given  by  cj  =  V(X+2p)/p. 

It  can  be  shown®^  that  the  displacement  vector,  u,  is  given  by  an  irrotational  potential, 
and  a  rotational  potential,  y.  For  the  case  of  an  infinitely  long  cylinder,  the  nonvanishing 
component  of  the  vector  v  is  y, ,  its  z-component.  Each  potential  satisfies  the  Helmholtz 
equation.  In  summary,  we  have: 

u  =  V(J>  +  Vxy  (2-8) 
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ar*  rdr  i^ee^ 


?li  l^y, 

dt^^r  dr  r^  d£^ 


v^v*=o 


(2-9) 

(2-10) 


where  kj = (©/cj  and  k, = (©/cj.  Again,  the  time  dependent  term  is  suppressed.  The 

radial  displacement,  u,.  and  the  normal  and  shear  stresses  within  the  target,  T„and  th» 
respectively,  are  given  by: 


L  a^  C f 

[d^  Tdr  r^ae^J 


(2-11) 

(2-12) 

(2-13) 


The  solution  of  the  Helmholtz  equations.  Equations  (2-9)  and  (2-10),  can  be  given  as: 


<|»(r,e)  =  po  i  iXgn Jn(kdr)cosne  (2 - 14) 

nsO 

V/r,0) =po  i  i"e„h„  Jn(l^r)smn0  (2  - 15) 

nsO 

The  part  of  the  solution  involving  the  Hankel  function,  (see  Equation  (2-3)),  must  be 
excluded  here  due  to  singularity  at  the  center  of  the  cylinder  at  r = 0 .  The  modal  constants  g„ 

and  h„,  along  with  b„  in  Equation  (2-3)  are  to  be  solved  by  an  alternative  set  of  boundary 
conditions,  namely,  at  r = a: 


Tjr=.=-Pi(a»0)  (2-16) 

u/a,  0) = u,(a,  0)  (2-1 7) 

Tj.=.  =  0  (2-18) 

Formally,  a  set  of  simultaneous  equations  involving  these  modal  constants  results: 

dnl>„+<l„g„+d,3h„  =  A,  (2- 19a) 

d2)bn+d22g„+d23h„  =  Aj  (2- 19b) 

d3.b„+d33g„+d33h„  =  0  (2-19C) 
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where 


d„=j[(k.a)'Hi»(k,a)] 

(2 -20  a) 

d,,  =  2(k.a)j:(kja)-Pn’-<k.a)“]  J.(k,a) 

(2 -20  b) 

d„  =  2n[-0£.a)l'.(lt,«)+JA'‘)l 

(2 -20  c) 

d„  =  -(l£,i)Hl''(k,a) 

(2-20d) 

d22  =  (k^a)j'„(kja) 

(2-20e) 

(2-20  f) 

d3,=0 

(2-20g) 

djj  =2n[J„(k„a)-(k„a)j’(k,a)] 

(2-20h) 

djj  =  2(k.a)  j;.(k,a)-  [2n^-<k,a)']  J„(k.a) 

(2-20i) 

A,  =  |[(k.a)'J.(k,a)] 

(2 -21  a) 

Aj  =  (k|  a)  J„(k,  a) 

(2 -21b) 

Thus  by  Crammer’s  rule,  b„is  given  by: 


(2-22) 


where  B„  and  D,  are  the  determinants: 
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d|2  di3 

B„  =  Aj  djj  djj 

0  dji  djj 

<^ii  d,2  d,j 

Da=  dj,  djj  djj 

0  d,2  djj  1 


(2 -23  a) 


(2 -23  b) 


In  the  previous  equations,  quantities  within  the  parentheses  are  the  arguments  of  the 
Bessel’s  functions.  Upon  expansion  of  the  determinants,  its  can  be  shown  that: 


J»(x)F,-xJ',(x) 

"  Hl'>(x)F.-xHl'y(x) 


(2-24) 


where  F.  is  given  by; 


F„=F.(x)=^x^ 


<122 

<*23 

<l22 

<*33 

<i.2 

<*.3 

<i22 

<*33 

(2-25) 


In  the  above  equations,  x  =  k,  a  and  Xt = k.a.  The  function  F„  is  a  measure  of  the  modal 
mechanical  admittance.  u,/t„,  of  the  elastic  cylinder.  This  function  contains  variables,  xx  and 
xa  (=  kja)  which  can  be  related  to  x  through  the  Lame’s  constants,  thus  in  effect,  F„  is  a 
function  of  x.  For  a  solid  with  no  damping,  F„  is  real.  Equation  (2-24)  reduces  to  Equation  (2-5) 
as  the  mechanical  admittance  of  the  target  becomes  nil  for  rigid  boundaries. 

In  the  far  field  where  k,  r  is  large,  the  Hankel  function  of  the  first  kind  can  be  given  by  its 
asymptotic  form: 


The  scattering  term  in  Equation  (2-3)  can  thus  be  given  by: 


(2-26) 


p*(k, r,0) = po 8„b„cosne  (2 - 27) 

For  subsequent  discussions,  it  is  convenient  to  introduce  a  modal  form  function  for  the 
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partial  waves,  namely: 


f.(x  ,0) = -j4=e„b«cosne 
VlTlX 

Equation  (2-27)  can  then  be  rewritten  as: 


4(x.e)= 


=  I  f.(x,e)  = 

asO 


2  * 

-77=  Z  e„b„cosne 

•ViJlx«»o 


(2-28) 


(2-29) 


where  fj(x,0)  is  the  normalized,  far-field  form  function  for  the  scattering  process. 

In  experimental  investigations,  the  scattered  power  commands  the  utmost  interest  and 
can  be  presented  by  the  bistatic  or  differential  scattering  crosssection,  o',  it  is  the  square  of 
the  modulus  of  the  normalized  form  function.  For  the  solid  cylinder  of  infinite  length,  this 
quantity  is  given  as: 


~-a'=IUx,er= 


S  f.(x,e) 

fisO 


2  « 

-r==  £  6,b,cosne 

■Vijrx»*o 


(2-30) 


For  backscattering  at  e  =  ji.  Equations  (2-28)  and  (2-29)  become: 


f„(x,n)=-=(-l)“e„b„ 

yinx 

f«(x  ,K)  =  I  f„(x  ,n) = -jl=(-l)Xb„ 
"=o  yin\ 

and  the  (normalized)  back-scattering  cross-section,  og,  is  given  by: 


na 


I  fn(x,x) 

ll«0 


IrtXU'O 


2:(-i)"£„b„ 


(2-31) 

(2-32) 


(2-33) 


Similar  equations  for  rigid  targets  can  likewise  be  derived.  Some  of  these  equations  for 
the  Infinite  cylinder  are  given  below  with  the  appropriate  superscript,  (r),  designation.  The 
modal  amplitudes  b„  have  already  been  given  in  Equation  (2-5).  Thus: 


(2-34) 

(2-35) 
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The  bistatic  scattering  and  the  backscattering  crosssections  are  respectively: 


(2-36) 

(2-37) 


The  clue  for  the  development  of  the  RST  was  prompted  by  the  realization  that  the 
scattering  crosssections  for  the  fluid  loaded,  elastic  target,  as  presented  by  the  magnitude  of 
the  far-field  form  functions,  are  very  similar  to  those  for  the  rigid  target  of  identical  geometry  with 
the  exceptions  of  frequency  locations  (k,a)  values  very  nearly  equal  to  the  natural  frequencies 

of  the  rigid  target.  The  similarity  (or  dissimilarity)  is  even  more  striking  when  individual  partial 
waves  for  the  two  classes  of  boundaries  are  compared  as  shown  in  Figure  2-2. 


FIGURE  2-2.  BACKSCA7TERED  SPECTRA  WITH  AND  WITHOUT  CONSIDERATION  OF 
RESONANCE^ 
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In  Figure  2-2(a),  amplitudes  of  several  partial  waves  from  (n  =  0  n  =  5)  for  a  rigid 
cylinder  are  plotted.  These  curves  are  given  by  Equation  (2-35)  and  are  a  set  of  smooth  curves 
as  compared  to,  in  Figure  2-2(b).  a  similar  set  of  more  spiky  ones  for  the  elastic  cylinder  from 
Equation  (2-31).  Ths  spikes  occur  at  frequency  locations  corresponding  to  the  free  vibration 
frequencies.  The  latter  are  the  real  parts  of  the  eigenfrequencies  associated  with  the 
determinant  D„  in  Equation  (2-23).  These  locations  are  marked  in  Figure  2-2(b).  in  parentheses, 
by  an  additional  index  /. 

Figures  2-3(a)  and  2-3(b)  depict  an  overlaid  comparison  of  the  two  sets  of  curves  taken 
from  the  previous  figure,  for  n  =  2  and  n  =  3,  with  their  respective  resonances  marked  by  the 
same  indices.  Thus  it  should  be  construed  that  In  Figure  2-3(a)  for  n  =  2  and  for  the  elastic 
cylinder,  there  are  resonances  (2,1),  (2,2),  (2,3),  etc.,  which  are  the  solutions  of  Dj  =  0  . 


FIGURE  2-3.  COMPARISON  OF  RIGID  AND  ELASTIC  RESPONSES"® 


Figure  2-4  gives  a  plot  of  the  difference  of  the  two  sets  of  curves  in  Figure  2-3(a)  or 
I  f2(rc,k,  a)  -  a)  |,  and  the  contention  that  the  scattered  pressure  for  a  fluid  loaded,  elastic 

target  can  be  separated  Into  a  rigid  background  and  a  resonance  component  (so  marked  by 
doubled  indices)  is  borne  out  in  an  even  dramatic  manner.  The  difference  is  due  to  the 
elasticity  or  penetrability  of  the  target  at  resonances  and  the  resulting  circumferential  waves  are 
referred  to  as  the  Rayleigh  surface  waves  for  /  =  1  and  as  whispering  gallery  waves  for  /  >  1 . 
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Figure  2-5(a)  is  an  overlaid  comparison  of  the  far  field  form  function  for  the  two  classes 
(rigid  and  elastic)  of  cylindrical  targets  as  a  function  of  the  angular  orientation.  Backscattering 
occurs  at  0  =  re.  The  form  functions  presented  are  approximations  containing  terms  of 
n = 0  -4  n = 22  at  k,  a  =  12.53  which  is  the  resonance  location  of  the  (3,2)  mode. 


FIGURE  2-5.  ASPECT  DEPENDENCE  OF  ELASTIC  AND  RIGID  RESPONSES 
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A  partial  expansion  of  the  form  function  for  the  elastic  cylinder,  taken  from  Equation  (2-30) 
is: 


|f.(I2.53,0)|  = 


|f.(12.53,0)l  = 
2 


22 


I  f„(12.53.0) 


+2 


Vl2.53n 
JjFj -12.53/2 


If,(12.53,0)_  _ 

n  =  0  I  |i>  =  0 

JoFo-12.53J’o  J,F, -12.53  j; 

— :  -  -  i-^cosQ 


H^'>Fo-12.53  H\‘>F,  -12.53  H 


(ly’ 


^  JjF, -12.53/3 

cos20  +  2-:;7:- — cos  30 


H<2''F2-12.53  -12.53  Hi'>' 

J22F22--I2.53/22 


■  +  2 


tf^>F22-12.53H^i^ 


cos220| 


(2-38) 


In  the  above  equation,  the  Bessel  and  Hankel  functions  and  their  derivatives  are 
computed  at  the  argument  of  k,  a  =  12.53 .  The  functions  Fo  through  F22  are  taken  from  Equation 

(2-25),  which  is  yet  another  set  of  Bessel  functions  of  various  orders  and  their  derivatives  at 
k,  a  =  12.53.  There  is  a  similar  series  for  the  rigid  cylindrical  target  which  is  the  same  as  that 

given  above  with  the  exception  that  all  the  admittance  functions  F„  =  0  for  (n = 0  -♦  22).  We  have; 


|f^\l2.53,0)|  =1 
2 


VII537C 


£  fi'’(12.53,0)  ~  Z  f^\l2.53,0) 

n=0  n=0 

j|;+2Acose  +  2icos2e 


+2j|jcos39+-.+2^cos22e| 


(2-39) 


As  seen  from  Figure  2-5(a),  there  are  obvious  differences  between  the  two  far-field  form 
functions  at  several  angular  orientations  including  the  backscattering  direction.  These 
differences  are  due  to  the  elasticity  of  the  target,  if  the  fourth  term  (n  =  3)  in  Equation  (2-39)  is 
replaced  by  a  similar  (fourth)  term  for  the  elastic  case  in  Equation  (2-38),  a  hybrid  equation 
designated  by  the  superscript  (h)  in  Equation  (2-39a),  results: 
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11^^12.53,0)1  = 


I  ^(12.53,0) 

usO 


22 


I  f2‘\l2.53.0) 

Q  zO 


Vl2.53re 


2  Jocose +  2icos2e 


^  J3F3-12.53J; 


■+2j^cos2201 


(2 -39  a) 


When  this  hybrid  equation  is  compared  to  that  of  the  rigid  case  in  Equation  (2-39)  as 
shown  in  Rgure  2-5(b),  the  differences  seen  in  Rgure  2-5(a)  more  or  less  disappear. 
Consideration  of  the  result  in  Rgure  2-5  and  that  in  Rgure  2-3  presented  previously  strongly 
leads  to  the  following  points: 

1 .  It  seems,  from  looking  at  Rgures  2-5(a)  and  2-5(b).  that  the  differences,  at 
resonance,  between  the  scattered  pressure  for  the  rigid  and  for  the  elastic  cases  are 
primarily  due  to  the  existence  of  a  single  elastic  term,  f3(l2J3,0).  pertinent  to  the  (3.2) 

resonance  and  the  (3,2)  resonance  alone.  This  implies  that  the  other  resonances  are  not 
at  all  excited  at  k, a  =  12.53. 

2.  It  also  implies  that,  term  for  term,  the  elastic  partial  wave  must  virtually  be  equal 
to  that  for  the  rigid  case  because  of  the  close  resemblance  of  the  elastic  and  the  hybrid 
curves  in  Rgure  2-5(b).  The  exception  is  the  f3(12.53,0)  term  mentioned  above. 

3.  In  Rgure  2-3,  the  differences  of  the  partial  waves  for  the  two  cases  were  shown 
for  the  backscattering  direction,  and  these  differences  are  reflected  in  the  (3,2) 
resonance  bistatic  form  function  at  0 =n  in  Rgure  2-5(a),  together  with  the  differences  at 
other  angular  orientations.  The  fact  that  a  single  elastic  resonance  term  f3(12.53,0) 

accounts  for  the  differences  not  only  for  0 = tc  but  for  all  other  angular  orientations,  see 
Rgure  2-5(b),  strongly  suggests  that  differences  in  partial  waves  are  not  limited  to  the 
backscattering  direction  but  extend  to  all  angular  orientations. 

From  the  above  discussion,  it  can  be  said  with  reasonable  conviction  that  the  scattered 
pressure  from  an  elastic  target  consists  of  a  rigid  background  due  to  specular  reflection  and 
diffraction  effects  and,  additionally,  a  resonance  component.  The  latter  is  due  to  the  elasticity 
or  penetrability  of  the  target  which,  by  virtue  of  admitting  energy,  gives  rise  to  the  many  modes 
of  free  vibration  within  the  target.  This  energy  admission  occurs  only  at,  and  very  near  one  of 
the  modal  resonances  of  the  target.  There  is  no  energy  admission/penetration  at  other 
frequencies,  and.  in  such  situations,  the  scattered  pressure  can  be  derived  as  if  the  target  is 
perfectly  rigid.  This  is  the  qualitative  statement  of  the  RST. 
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The  quantitative  statement  of  the  RSI  involves  the  explicit  subtraction  of  the  rigid 
background  from  the  elastic  form  function.  Defining  an  elastic  scattering  function  S„(x)  such 

that: 


S,(x)-l=2b„  =  -2 


J,(x)F.-xJ',(x) 

Hl«(x)F„-xH<>’^(x) 


and  further  defining  an  elastic  scattering  phase  shift  ii„(x)  such  that: 


taiiTi„(x)  = 


J,(x)F,-xj;(x) 

N„(x)F„-xK(x) 


(2-40) 


(2-41) 


where  N|.(x)  is  the  Bessel  function  of  the  second  kind  and  relates  to  the  Hankel  function  by 
H«^(x) = J.(x) + iN,(x).  The  conjugate  of  is  H?>(x) = J.(x)  -  iN,(x). 

It  can  be  shown  by  taking  the  reciprocal  of  b„(x)  from  Equation  (2-24)  and  simplifying  that: 


b„  =  (2  -  42) 

and  by  substituting  b„(x)  back  into  Equation  (2-40),  we  have: 


S„=  1 +2b„=  1 +2ismii„e  =e 

(2-43) 

A  set  of  equations  similarly  exists  for  the  rigid  case: 

s;"(x)-l=2b«  =  -2^(x) 

(2-44) 

.  (rV  -v  ^n(x) 

(2-45) 

b?  =  ismTi®J‘’^^ 

(2-46) 

(2-47) 

Also  from  Equation  (2-44), 

^  to  to-iKCx) 

•  Hi'^Cx)  »s\x)  ~  H<'>'(x) 

(2-48) 
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We  now  want  to  split  the  elastic  scattering  function,  Sb(x)  into  a  rigid  portion  and  a 
resonance  portion: 


Hg>(x)F,(x)-xHg^(x) 

H<‘>(x)F.(x)-xH<‘’^(x) 


(2-49) 


It  should  now  be  recalled  that  Fa(x)  is  a  measure  of  the  admittance  of  the  target,  hence  its 
reciprocal  can  be  viewed  as  the  modal  mechanical  impedance  of  the  target.  The  conjugate 
terms  H®/x  H®,  (i  =  1 , 2)  can  be  viewed  as  the  modal  acoustical  impedance  of  the  fluid  given 

by: 


Hi*\x) 

xH?>'(x) 


=  A„(x)+iK„(x); 


"^^=A,(x)-iK.(x) 


;HO''(x) 


where: 


1J„(x)J;(x)+N„(x)n;(x) 

A.(x)=-- 


K„(x)  =  - 


*  [/.(xrf+lKWf 

2  I 


"’‘‘[I’.Wf+tK.Cxrf 


The  quantity  K,(x)  in  Equation  (2-52)  is  negative. 

Let  a  resonance  frequency,  x„/,  be  defined  such  that: 


lxHf(x.,) 


=  ^,(x„;) 


(2-50) 


(2-51) 

(2-52) 


(2-53) 


The  above  characteristic  equation  is  equivalent  to  the  statement  attributed  to  Junger  and 
FelP  that  resonances  should  appear  when  the  sum  of  the  mechanical  and  the  radiation 
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impedances  goes  to  zero.  In  Equation  (2-53),  the  resonance  frequency  is  designated  by 
indices  n  and  /  witti  the  index  n  pertinent  to  each  of  the  partial  waves.  The  additional  index  / 
indicates  the  multiplicity  of  the  solution,  namely,  for  each  partial  wave,  there  exist  the 
fundamental  resonance  (n,  1)  and  harmonics  (n,/),  theoretically  infinite  in  number.  Their 
frequencies  are  very  nearly  equal,  but  not  identical,  to  the  (complex)  eigenfrequencies  obtained 
from  the  Equation  D„(x)  =  0. 

We  now  adopt  the  linear-approximation  method  of  the  nuclear  resonance  theory®®  by  expanding 
the  impedance  term  in  the  vicinity  of  one  of  the  resonance  frequencies,  x„„  in  a  Taylor  series: 


f;;'(x)  -  f;'(x„,) = f;'(x)  -  A„(x„;)~p„,(x  -x„,) 


(2-54) 


By  substituting  the  results  in  Equations  (2-48).  (2-54),  and  (2-50)  into  Equation  (2-49)  to  obtain: 


S.(x>-S'V) 


x-x„,+5ir„; 


(2-55) 


The  quantities  and  r„,  are  undefined  constants  and  are  related  by: 


2K„/ 

(2-56) 

Pn/ 

The  quantity  r„,  in  Equation  (2-56)  is  positive  and  can  be  viewed  as  the  width  of  the  (n,/) 
resonance.  Thus: 


s.(x)-i=.sr(x)— 

x-x„,+-iT„, 


=  S^) 


1  — 


ir. 


n/ 


1  •  1 


x-x„,+jir„, 


-1 


=  (^'’(x)-l]-S”(x) 


ir, 


n/ 


x-x„,+5ir„, 


(2-57) 


and  by  means  of  Equation  (2-47), 
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S,(x)- 1 


x-x„,+i|T„, 


(2-58) 


The  summation  over  /  =  1  ->  oo  is  added  in  the  last  expression  due  to  the  multiplicity  in  the 
solution  of  Equation  (2-53)  as  explained  previously.  This  summation  generally  does  not  refer  to 
the  summing  of  infinite  terms  but  rather  it  denotes  the  activation,  one  term  at  a  time,  as  one  of 
the  resonances  x„  is  approaching. 


The  substitution  into  Equatior'  (2-28)  yields  the  final  form: 


2  1 

^(x  ,6) = -j==£„b„cosne = -===e„[S„(x)  - 1]  cosnG 
VlJtX  \17CX 


-jr=£„2  sinti® +  I 
yinx 


'='x„,-x--ir„, 


cosnG 


(2-59) 


This  equation  explicitly  shows  that  the  form  function,  and  the  scattered  pressure  thereof, 
given  by  Equation  (2-30),  consists  of  a  rigid  background  and  a  resonance  portion  which  comes 
into  effect  as  the  resonances  of  the  partial  waves  approach.  It  constitutes  the  quantitative 
statement  of  the  RST. 


The  mathematical  treatment  for  other  cases  is  similar  to  that  ior  the  infinite,  solid  cylinder. 
For  a  hollow  shell,  two  unknown  constants  associated  with  the  Hankel  function  are  to  be 
included  In  Equations  (2-14)  and  (2-15).  An  equation  similar  to  Equation  (2-14)  to  describe  the 
pressure  field  in  the  enclosed  fluid  exists  with  an  extra,  unknown  constant.  This  brings  the 
unknown  constants  to  a  total  of  six  and  results  in  6  by  6  determinants  in  B„(x)  and  D„(x)  in 

Equation  (2-22).  An  additional  layer  of  shell  will  increase  the  order  of  the  matrices  by  4  to  10. 
For  spherical  cases,  the  elements  of  the  matrices  will  have  become  spherical  Bessel  functions 
and  the  angular  variable  of  the  scattering  is  dictated  by  the  Legendre  polynomials  with  (cos 6) 
as  their  argument. 
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CHAPTER  3 

EXPERIMENTAL  FACILITY  AND  ARRANGEMENT 


The  experimentation  to  verify  the  validity  of  the  RST  and  to  develop  the  methodology  for 
underwater  target  classification  by  means  of  active  sonar  was  conducted  in  the  Naval  Surface 
Warfare  Center  Coastal  System  Station  (NSWCCSS)  Hydroacoustic  Measurements  Facility  at 
White  Oak,  Maryland.  The  facility  is  one  of  the  two  mayor  testing  facilities^^-^  operated  and 
managed  by  the  Mine  Systems  Branch  (Code  2620)  of  the  Mine  Systems  Division  (Code  260) 
which  is  a  division  of  the  Costal  Warfare  Systems  Department  (Code  20). 

The  White  Oak  facility,  located  in  Building  217,  is  equipped  with  a  redwood  test  tank  30’ 
in  diameter  and  20’  in  depth.  This  tank  is  placed  underground  within  a  36’  diameter,  steel, 
earth-retaining  drum  which  sits  on  a  resilient  support  on  a  2’  reinforced  concrete  foundation. 
Solid  earth  around  the  drum  and  the  concrete  foundation  effectively  isolate  the  tank  on  all  sides 
except  the  top  side  from  building,  roadway  and  other  operating  equipment  vibration  and  noise. 
The  arrangement  also  provides  for  electrical  shielding.  Built  over  the  tank  is  the  entire 
church-like  building  which  is  constructed  witti  laminated  beams  with  a  very  high-peaked  roof  to 
facilitate  a  precision  hoisting  or  lifting  carriage.  The  main  floor  or  the  test  bay,  which  is  4’  above 
the  top  of  the  tank,  is  made  of  12",  post-stressed,  reinforced  concrete  with  a  52"  wide  by  46’  slot 
running  over  the  diameter  of  the  tank.  Heavy  covers  over  the  slot  are  provided  to  seal  off  the 
tank  from  the  outside  resulting  in  an  isothermgd  testing  environment  and  further  mitigating 
airborne  noise  intrusion.  The  tank  is  filled  with  continuously  filtered  fresh  water.  The  filtering 
process  not  only  minimizes  unwanted  scattering  from  floating  or  suspending  particles  but  also 
affords  a  refraction-free  test  environment  by  eliminating  temperature  gradients  in  the  water.  The 
parking  lot  is  at  the  same  level  as  the  test  bay,  and  equipment  weighing  as  much  as  8,000 
pounds  can  be  placed  on  the  floor  and  moved  around  using  a  folk  lift.  A  transducer  supporting 
and  positioning  structure  is  isolated  mechanically  from  the  building  and  the  tank.  It  consists  of 
a  system  of  interconnecting  vertical  columns  resting  on  long,  horizontal  "grade"  beams 
underground.  On  the  top  of  the  structure  are  65’  long  horizontal  guide-rails  (crane  rail  beams) 
which  support  and  guide  a  lifting  carriage  with  a  two-ton  capacity.  This  lifting  carriage  is 
capable  of  handling  all  moving  requirements  wth  speed  and  precision  including  horizontal 
movements  along  the  tank  diameter  and  vertical  lifting  and  positioning  of  loads  or  transducers. 
A  cross  section  of  the  building  and  the  tanks  is  given  in  Figure  3-1 . 

At  the  test  bay  level  are  two  "I"  beams  similar  to  the  crane  rail  beams  which  are  also 
supported  independently  and  acoustically  isolated  from  the  vertical  columns  and  the  floor. 
These  beams  support  the  slot  covers  and  two  additional  transducer  supporting  and  positioning 
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platforms  which  can  be  moved  by  sliding  along  the  slot  on  circular  guides  or  by  lifting  and 
lowering  to  position  by  the  lifting  carriage.  A  standard  transducer  platform  has  fixtures  which 
can  place  a  transducer  or  projector  to  any  depth  within  the  tank.  Precise  positioning  along  the 
slot  is  accomplished  by  a  motorized  drive.  The  other  platform  which  has  a  single  axle  that 
rotates  above  the  tank  and  parallel  to  the  water  surface.  Attached  to  the  end  of  the  axle  is  a 
horizontal  positioner  for  the  suspension  of  hydrophones  or  loads  at  a  desired  test  level  in  water. 
The  arrangement  is  suitable  for  rotating  a  target  in  a  horizontal  plane  about  a  designated  axis 
as  Is  used  to  determine  the  directivity  patterns  of  nonsymmetrical  line  hydrophones.  The 
rotating  axis  is  made  of  a  hollow  tube  which  also  allows  a  target  to  be  suspended  axially,  but  at 
a  desired  angle,  0,  with  respect  to  the  incoming  wave,  while  a  hydrophone  riding  on  the 
positioner  rotates  around  the  target  at  a  distance,  R,  from  the  target  as  depicted  in  Figure  2-1. 
Both  platforms  slide  on  two  round  tracks,  each  46’  long,  which  are  level  and  at  a  constant 
separation  for  precise  alignment. 


STEEL  TANK 


■88'  8' 


FIGURE  3-1.  WHITE  OAK  SITE  OF  HYDROACOUSTIC  MEASUREMENTS  FACILITY 


The  heart  of  the  signal  generation,  conditioning,  and  recording  rests  within  an  electronic 
console  which  is  housed  in  a  separate  air-conditioned  room  on  the  main  floor.  The  unit  has  a 
transmitting  section  which  supplies  power  to  the  projector  and  a  corresponding  receiving 
section  which  accepts  and  modifies  transduced  signal  from  the  hydrophone.  Both  sections  are 
controlled  by  a  pulse  time  generator  which  sets  the  transmitting  pulse  width,  the  receiving  pulse 
width  and  its  time  delay  with  respect  to  the  transmitted  pulse.  Operational  frequency  varies 
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continuously  from  50  Hz  to  2  MHz  with  a  50-clB  dynamic  range.  The  synthesizer  generates  CW 
signal  adjustable  from  1  to  99  seconds  over  a  pre-determined  frequency  range  in  linear  or 
logarithmic  scale.  The  CW  signal  is  broken  up  by  the  opening  and  closing  actions  of  the 
transmitter  signal  gate  into  pulses  which  eventually  are  emitted  through  the  projector  as  a  train 
of  pressure  pulses  in  water.  The  (transmitted)  pulse  width.  T„  and  the  rapidity  of  the  actions  or 

repetition  rate  are  set  by  the  pulse  timing  generator.  These  pressure  pulses  are  sampled  by  the 
hydrophone  at  a  time  delay,  T^,  with  respect  to  the  opening  of  the  transmitter  signal  gate  and  at 

a  sampling  (received)  pulse  width  T,.  Both  quantities  are  also  set  by  the  pulse  timing  generator. 

For  each  pulse,  no  signal  will  be  detected  by  a  transducer  beyond  (Tj+T,).  Recording 

instruments  within  the  electronic  console  include  a  rectangular  recorder  with  a  linear  or 
logarithmic  paper  drive  and  a  polar  recorder  which  are  pegged  automatically  to  the  various 
functions  and  controls  of  the  synthesizer,  the  pulse  timing  generator  and  the  rotating 
mechanism.  Analog  signal  at  various  points  in  the  console  is  monitored  by  dual-beam  storage 
oscilloscopes  with  delay  circuitry.  Pictures  of  the  pulse  at  various  stages  of  resonance  can  be 
recorded  in  photographs.  A  block  diagram  of  the  major  components  of  the  electronic  console 
is  given  in  Figure  3-2. 

RST  experiments  have  been  done  exclusively  in  the  White  Oak  Site  in  the  large  tank.  The 
other  major  testing  facility  at  Brighton  Dam,  about  14  miles  north  of  White  Oak,  Maryland,  offers 
a  comparable  facility  in  a  barge  on  an  800-acre  reservoir.  Average  test  depth  is  25’.  Future 
experiments,  especially  on  CW  mode,  can  be  performed  in  the  Brighton  Site. 

In  summary,  the  electronic  console,  the  moving  carriage  and  the  two  platforms  allow  for  a 
variety  of  projector-target-hydrophone  (P-T -H)  orientations  and  many  types  of  experimentation 
for  the  RST  investigation.  Some  of  the  notable  possibilities  are  given  as  follows. 

Frequency  response  over  a  predetermined  frequency  range  in  linear  or  log  scale  is  very 
easily  obtained.  Figure  3-3  depicts  the  data  from  this  type  of  tests. 

For  this  experiment,  the  hydrophone  is  typically  positioned  between  a  projector  and  its 
target  in  a  line,  and  in  the  middle  of  the  tank  at  mid-depth.  Depending  on  the  time  delay,  the 
recorded  signal  can  be  the  direct  pulse  which  insonifies  the  target,  as  shown  in  the  top  trace, 
the  backscattering  or  the  monostatic  cross  section  of  the  target  as  the  middle  trace,  or  at  a 
further  time  delay,  the  reradiated  of  the  target  with  the  background  suppressed  as  the  bottom 
curve.  The  latter  two  curves  are  expected  to  yield  the  various  modal  resonances  as  speculated 
upon  in  the  introductory  section. 

The  previous  description  of  the  data  collection  system  reflects  what  was  used  to  collect 
the  data  in  this  report.  It  has  since  been  upgraded  to  a  computer  controlled  data  collection 
system. 
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FIGURE  3-3.  SPECTROGRAM  DATA 
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Monostatic  broadside  spectrograms  can  be  measured  using  the  test  geometries  shown  in 
Figure  3-4. 


If  the  frequency  is  slowly  varied  by  taking  smaller  frequency  steps  over  less  of  a  span  a 
"Zoom"  effect  can  be  achieved  and  refined  features  near  and  at  a  modal  resonance  such  as  the 
bandwidth  of  the  resonance  peak,  and  its  shape  {the  rapidity  or  slope  at  which  a  resonance  is 
arrived  or  receded)  can  be  obtained. 

Bistatic  frequency  response  can  be  obtained  with  a  P-T-H  arrangement  shown  in  Figure 
3-5  The  target  can  either  be  in  a  vertical  or  a  horizontal  orientation.  When  the  target  is  in  the 
vertical  orientation,  the  insonification  on  the  target  remains  normal,  while  oblique  incidence  on 
the  target  will  result  when  the  target  is  in  a  horizontal  orientation. 

The  aspect  dependence  of  the  bistatic  scattering  cross  section  can  be  obtained  with  the 
target  suspended  freely  at  the  axis  of  the  rotator  and  the  hydrophone  attached  to  the  rotating 
positioner.  The  resulting  directivity  pattern  can  be  at  any  frequency  including  resonance 
frequencies.  The  geometry  for  this  test  is  shown  In  Figure  3-6. 


FIGURE  3-5.  BISTATIC  TEST  CONDITION  FOR  SPECTROGRAMS 
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FIGURE  3-6.  TEST  CONDITIONS  FOR  ASPECT  DEPENDENCE  MEASUREMENTS 

As  shown  the  target  can  be  either  in  the  vertical  or  horizontal  orientation.  The  hydrophone 
or  the  target  can  be  rotated  while  the  other  is  held  stationary  at  any  arbitrary  angle  or  both  can 
be  rotated  in  synchronous  fashion. 

One  of  the  most  important  features  of  any  acoustic  measurements  system  is  the  ability  to 
time  gate  the  hydrophone  signal  to  measure  the  desired  signal  and  to  reject  unwanted  noise. 

In  a  closed  environment,  one  source  of  noise  is  the  unwanted  reflections.  The  best  way  to  set 
these  delay  times  and  to  find  unexpected  echos  is  to  display  the  pulses  on  an  oscilloscope.  It 
is  possible,  with  the  proper  settings,  to  view  the  direct,  reflected  and  the  re-radiated  pulses 
either  individually  or  all  at  one  time.  Slow  frequency  sweeps  can  then  be  made  to  insure  that 
the  test  geometries  and  delay  selections  are  correct.  Then  the  measurements  become  routine. 
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Figure  3-3  has  three  data  curves  on  it  which  were  collected  by  setting  the  receive  gate  at 
three  different  delay  times  referenced  to  the  transmission  of  the  pulse.  The  upper  curve  is  the 
direct  pulse  which  provides  a  means  of  determining  the  intensity  of  the  signal  reaching  the 
target  or  being  transmitted  from  the  projector,  it  was  collected  by  setting  the  receive  gate  at 
the  delay  time  (plus  a  little  extra)  ttiat  it  takes  the  signal  to  go  from  the  projector  to  the 
hydrophone.  The  middle  curve  is  a  measure  of  the  intensity  of  the  signal  reflected  from  the 
target  and  is  dependent  on  the  differential  scattering  cross-section  pattern,  as  described  by 
Equation  (2-30).  For  this  measurement  the  receive  gate  must  be  set  toward  the  end  of  the 
reflected  pulse.  The  bottom  curve  is  a  measurement  of  the  intensity  of  the  re-radiated  energy. 
For  this  measurement  the  receive  gate  must  be  set  toward  the  beginning  of  the  re-radiated 
pulse.  This  is  easy  to  see  in  Figure  3-7. 


FIGURE  3-7.  DEPICTION  OF  THE  TIME  DOMAIN  SIGNAL 

in  the  real  world  the  transmitted  pulses  are  not  square  since  it  takes  a  cycle  or  so  for  the 
transducer  to  attain  equilibrium.  There  is  a  corresponding  number  of  cycles  at  the  end  before 
the  transducer  comes  to  rest.  This  means  that  there  should  always  be  a  cycle  or  two  between 
the  leading  edge  or  the  trailing  edge  of  the  pulse  and  the  sampling  gate  in  order  to  eliminate 
these  unstable  portions. 

The  arrangements  discussed  up  to  this  point  pertain  more  or  less  to  an  Infinite  cylinder 
suspended  freely  in  the  water.  Of  course  the  cylinders  we  measured  had  finite  length  but  since 
the  tests  were  conducted  using  signals  in  the  hundreds  of  kilohertz  region  our  projectors  had 
finite  beam  width  and  as  long  as  only  a  portion  of  the  cylinder  was  insonified,  it  reacts  like  an 
infinite  cylinder.  Rotation  of  the  cylinder  by  itself  creates  no  new  situation  due  to  the  symmetry 
of  the  target.  In  the  vertical  orientation,  and  the  impingement  was  always  normal  to  the  target. 
The  same  applies  for  a  target  which  is  spherical.  However,  the  situation  changes  for  other 
targets  whose  surface  is  not  symmetrical  to  the  (vertical)  axis  of  suspension.  An  example  is  a 
finite  cylinder  suspended  with  its  axis  in  a  horizontal  orientation.  Another  example  is  a 
horizontally  oriented  spheroid  with  its  circular  cross  section  axis  perpendicular  to  the  free 
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surface.  With  the  present  test  facility,  the  rotation  of  the  target  itself  and  the  independent 
rotating  of  the  hydrophone  on  the  horizontal  positioner  create  many  additional  testing 
possibilities. 
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CHAPTER  4 

RST  DEPLOYMENT  METHODOLOGY 


Before  a  methodology  for  the  detection  and  classification  of  underwater  targets  can  be 
finalized  for  deployment  at  sea.  it  is  necessary  to  develop,  within  the  context  of  the  available 
facility  and  its  capability,  experiments  which  \Mli  verify  the  predictions  of  classic  resonance 
theory.  As  mentioned  in  the  introductory  section,  successful  validation  would  automatically 
provide  the  sought-for  methodology  for  the  active  classification  of  underwater  targets.  At  the 
least,  the  experience  and  expertise  gained  by  the  venture  will  be  invaluable  assets  for  future 
methodology  development  and  improvement 

in  the  planning  of  the  experiments,  the  prevention  of  wave  interference  commanded  the 
most  attention.  This  consideration  led  to  the  use  of  the  pulse  mode  for  the  experimentation; 
and.  in  turn,  to  the  many  decisions  pertaining  to  the  physical  arrangements  of  the  testing 
instruments,  their  sizes,  and  to  the  ways  and  timing  the  signals  were  sampled,  etc.  These 
considerations  are  discussed  in  the  following  paragraphs. 

One  of  the  first  considerations  to  address  was  the  placement  of  the  projector,  target,  and 
the  hydrophone.  Recalling  Figure  3-5  and  the  associated  discussion  the  bulk  of  the 
consideration  was  to  eliminate  the  unwanted  reflections.  The  target  should  preferably  be 
placed  in  the  central  region  of  the  tank  which  is  approximately  15’  from  the  tank  walls  and  at  a 
mid-depth  of  10’  from  the  free  surface.  The  projector  radiates  in  a  rather  narrow  beamwidth, 
which  is  inversely  proportios  al  to  the  driving  frequency.  It  should  be  placed  far  from  tiie  target 
to  assure  that  the  rays  impinging  on  the  target  are  reasonably  parallel  and  even  in  intensity  to 
simulate  a  plane  wave.  The  hydrophone  which  samples  the  reflection  should  also  be  placed 
reasonably  far  from  the  target  to  achieve  as  good  a  far-field  condition  as  possible.  A  good 
compromise  led  to  the  placement  of  the  target  at  the  center  of  the  tank  with  the  projector  10’ 
from  it.  The  hydrophone  was  placed  3.5’  from  the  target.  This  P-T-H  arrangement  was  used  for 
the  majority  of  the  experiments. 

The  next  consideration  was  the  selection  of  the  size  of  the  target,  hydrophone  and 
projector.  The  choice  of  the  target  size  is  contingent  upon  two  major  factors.  First,  it  must  be 
small  compared  to  the  physical  arrangement  described  in  the  last  paragraph  to  satisfy  (a)  the 
plane  wave  incidence  condition  which  allows  an  even  insonification  on  the  target,  and  (b)  the 
far  field  condition  with  respect  to  the  hydrophone  sampling  for  scattering  or  re-radiation. 
Secondly,  the  size  must  be  small  In  comparison  with  the  wavelength  over  the  frequency  (ka) 
range  of  the  theoretical  predictions.  For  most  of  the  experiments,  cylindrical  rods  from  0.75"  to 
1.87"  in  diameter  and  6*  or  longer  in  length  were  used.  Those  diameters  yielded  ka  values 
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ranging  from  4  (0.75"  rod  at  100  kHz)  to  50  (1 .87"  rod  at  50  kHz).  A  length  of  6"  was  long 
enough  to  minimize  end  effects  and  preserve  the  two-dimensional  scattering  condition  of  an 
infinite  cylinder. 

In  addition  to  the  cylindrical  rods,  other  targets  were  chosen  for  experimentation, 
including,  hollow  cylindrical  shells  of  various  sizes  (0.75"  to  6.0")  with  ID/OD  ratio  of  2/3  and 
ID/OD  ratio  very  nearly  equal  to  1.0  (very  thin-walled  cylindrical  shells);  an  ellipsoid  3.5*  by  0.7* 
measured  in  the  major  and  the  minor  axes;  a  brass,  hollow  spheroid  approximately  3.5*  by  7"  in 
outer  size;  and  a  4.5*  diameter,  steel  pressure  tank.  The  hollow  cylindrical  shells  of  aspect  ratio 
of  2/3  was  used  to  verify  previous  RSI  predictions  and  other  targets  were  included  to  add 
variety  to  the  experimentation.  These  models  represented  a  first  step  toward  more  practical 
targets.  A  photograph  of  several  of  the  targets  is  given  in  Figure  4-1 .  The  targets  were  made  of 
aluminum  unless  otherwise  stated.  The  material  had  not  been  mentioned  prominently  because 
it  played  very  little  role  In  the  experimental  verification  of  the  theoretical  prediction  since  the 
original  RST  formulation  was  developed  using  a  ratio  wave  number  in  the  material  to  the  wave 
number  of  the  medium. 

In  the  choice  of  a  hydrophone  for  the  experiments,  size  was  again  an  important 
consideration.  A  Celesco  LC-10  '  .yurophone  was  selected  mainly  because  its  diameter  of 
0.250"  offered  very  little  blockage  of  the  incident  wave  to  the  target  when  placed  in  the 
back-scattered  position.  The  loss  of  transducer  sensitivity  in  the  choice  of  a  small  transducer 
was  partially  compensated  for  the  gain  in  the  applicable  frequency  range.  It  is  also  advisable 
to  operate  below  the  resonance  frequency  of  the  hydrophone,  which  for  the  LC-10  Is  225  kHz. 

Size  was  a  less  critical  factor  in  the  choice  of  a  projector.  A  USRD  E-27  transducer 
having  an  effective  diameter  of  0.4"  was  chosen  because  its  operating  range  extended  from  50 
to  700  kHz. 

The  third  consideration  was  the  transmitting  and  receiving  gate  or  pulse,  width.  As 
discussed  in  the  experimental  testing  facility  section  and  graphically  illustrated  in  Figure  3-2,  a 
transmitted  pulse  starts  when  the  transmitter  signal  gate  of  the  pulse  timing  generator  is 
activated  (goes  to  a  logical  "1"  and  ends  when  the  same  gate  is  deactivated  (goes  to  "0").  The 
time  lapse,  T, ,  is  the  transmitted  pulse  width.  Similarly,  after  a  chosen  time  delay  with  respect 

to  the  activation  of  the  transmitted  signal  gate,  ,  the  received  signal  gate  is  activated  and 

deactivated  over  a  time  interval,  T, ,  during  which  a  portion  of  the  transmitted  or  reflected  signal 

as  received  by  the  hydrophone  is  sampled.  This  sequence  of  actions  is  repeated  many  times 
per  second  and  is  referred  to  as  the  pulse  repetition  rate.  Thus  a  train  of  pulses  can  be 
identified  for  each  transmitted  pulse,  and  the  time  between  these  pulses  is  the  difference  In 
time  of  flight.  The  difference  in  flight  path  is  given  by  the  product  of  the  difference  in  flight  time 
and  the  speed  of  propagation.  The  position  of  the  transducers  and  the  target  must  be 
optimized  to  provide  minimal  interference  at  the  time  and  position  of  the  measurements. 
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This  type  test  requires  large  source  levels  ar^d  thus  a  very  large  driving  voltage.  Since 
the  electronics  gating  circuit  does  not  switch  this  signal  to  zero  volt,  there  can  still  be  a  signal 
being  produced  by  the  projector  and  received  at  the  hydrophone  that  can  have  similar 
amplitude  to  the  signal  being  measured.  This  residual,  plus  the  fact  that  the  projector  is  in 
vigorous  motion  when  the  pulse  goes  to  “0"  causes  a  exponentially  decaying  tail  on  the  training 
edge  of  the  pulses.  Sometimes  these  tails  cannot  be  avoided  and  also  add  to  the  sampled 
signal.  Rgure  4-2  is  a  depiction  of  the  pulse  as  it  appears  in  the  water.  In  later  figures  the 
delay  time  to  the  leading  edge  of  the  pulse,  in  microseconds,  is  placed  on  the  pulse  depiction 
where  the  "XXXX*  is  in  this  figure. 


Rgure  4-2.  depiction  of  acoustic  pulse 


In  Figure  4-3,  the  locations  of  the  tank  walls,  back  (BW)  and  front  (FW),  and  the  projector 
(P),  hydrophone  (H),  target  (T),  and  traveling  pulse  are  shown  in  a  linear  scale.  The  transmitted 
pulse  of  constant  amplitude  is  depicted  by  its  rectangular  envelope  and,  based  on  a  sound 
speed  of  4.875  ft/ms,  has  a  length  of  4.875  ft.  The  ends  of  the  envelope  depict  the  leading  and 
the  trailing  cycles  of  the  wave  train.  The  length  of  the  "tail",  or  residue  is  taken  as  1  ft. 


FIGURE  4-3.  SAMPLING  OF  DIRECT  PULSE 


A  serie.*?  of  four  pictures  showing  (i)  the  leading  edge  of  the  transmitted  pulse  sinusoid 
reaching  the  hydrophone,  (ii)  some  designated  time  later  when  the  wave  train  is  first  sampled, 
(iii)  the  trailing  edge  leaving  the  hydrophone,  and  (iv)  a  very  short  time  later  when  the  residue  is 
sampled.  The  numbers  on  the  pulse  are  associated  with  the  delay  time  in  microseconds  as  if 
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the  picture  was  frozen  at  that  instance.  The  arrows  indicate  the  direction  of  propagation.  At 
this  early  moment,  no  reflection  exists  and  the  received  wave  train  (shown  darkened)  sampled 
by  the  hydrophone  in  case  (ii)  corresponds  to  the  direct  signal  and  is  a  measure  of  the  intensity 
of  the  incident  signal. 


In  Figure  4-4,  the  incident  wave  has  reached  the  target  and  reflections  are  shown  (by 
solid  rectangles).  Recalling  the  discussion  associated  with  Figure  3-7,  the  reflected  pulse  is 
sampled  near  the  trailing  edge  to  measure  the  maximum  effect  of  the  resonance  as  depicted  in 
case  (ii)  of  Figure  4-4.  The  amplitude  of  the  reflection  is  generally  much  smaller,  by  some  20 
db,  and  is  therefore  shown  pictorially  by  a  thinner  rectangular  envelope.  At  a  resonance,  the 
amplitude  will  even  be  smaller  due  to  energy  penetration  into  the  target  as  discussed  in 
Chapter  3.  Of  course,  the  sampling  location  within  the  envelope  is  Immaterial  for  the 
nonresonance  cases  since  the  envelope  remains  rectangular.  The  reduction  in  target  strength 
at  the  resonances  allows  one  to  determine  the  signature  of  a  particular  underwater  target  as 
shown  in  Rgure  3-3. 

The  result  of  sampling  in  case  (ii),  whether  away  from  resonance  or  at  a  resonance  state, 
is  due  to  the  sum  of  infinite  partial  waves.  Equation  (2-33)  from  the  theory  section  is  repeated 
here  for  clarity. 
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(2-33) 


The  incident  wave  scatters  all  modes;  and  as  long  as  there  is  wave  incidence  on  the  target, 
there  always  exist  scattered  waves  of  various  frequencies.  However,  if  sampling  is  slightly 
delayed  to  allow  the  trailing  cycles  of  the  reflected  wave  train  to  be  sampled  by  the  hydrophone 
as  shown  in  case  (iv),  the  sampling  result  will  be  that  due  to  the  reflection  from  the  residue 
transient  which  is  generally  of  negligible  magnitude.  The  exception  occurs  when  the  target  is  at 
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a  resonance  state.  At  resonance,  the  target  reradiates  its  stored  energy  at  a  single  frequency 
because,  as  predicted  by  the  resonant  theory,  only  the  energy  associated  with  the  resonance 
frequency  would  have  penetrated  into  the  target.  The  energy  reradiation  at  the  end  of  wave 
incidence  at  resonance  is  at  a  rate  controlled  by  the  damping  of  the  target  material  and  the 
radiation.  The  energy  reradiation  has  a  magnitude  much  larger  than  the  reflection  of  the 
residue  in  the  "tail":  and  if  sampling  is  made  in  this  manner  over  a  series  of  frequencies,  the 
result  is  called  the  spectrogram  of  the  reradiation  or  the  free  ringing  curve  with  peaks  at  the 
various  resonance  frequencies.  This  is  an  alternative  method  to  obtain  the  signature  of  the 
underwater  target. 

Figure  4-5  depicts  the  back  wall  reflection  from  the  incident  wave  coming  out  from  the 
back  or  baffled  side  of  the  projector.  The  amplitude  of  this  reflection  is  too  small  to  warrant 
serious  concern  pertaining  to  wave  interference. 


Figure  4-6  depicts  the  front  wall  reflection  of  the  forward,  incident  wave.  This  pulse 
passes  the  hydrophone  as  depicted  in  Figure  4-3,  proceeds  to  the  front  wall  of  the  tank  and 
comes  back  to  the  hydrophone  after  travelling  through  a  distance  of  43.5’.  The  amplitude  of 
this  reflected  pulse  is  comparable  to  the  first  echo  from  the  target  in  Figure  4-4,  but  because  of 
an  additional  travel  of  30’,  there  is  a  considerable  time  delay  in  arrival  at  the  hydrophone  and  it 
can  be  eliminated  by  time  gating. 

There  is  at  least  one  more  reflection  worthy  of  consideration  which  is  not  depicted  in 
these  figures.  That  is  the  reflection  from  the  free  surface  of  the  tank.  The  pulse  would  have 
travelled  a  distance  of  slightly  more  than  20’:  10’  up  from  the  projector  to  the  free  surface  and 
then  10’  down  to  reach  the  hydrophone.  The  corresponding  time  for  the  pulse  to  pass  the 
hydrophone  covers  approximately  a  period  from  4,300  to  5,300  ns.  The  amplitude  of  this 
reflection  is  variable  depending  on  the  directivity  of  the  projector  at  the  various  incidence 
frequencies,  and  it  may  be  comparable  to  that  of  the  backsc  ^ttered  reflection.  However,  it  is 
important  to  include  this  reflection  in  the  study  of  wave  intenerence  because  of  its  effect  on  the 
backscattered  reflection  and  the  ensuing  reradiated.  This  limits  the  delay  in  sampling  the 
reradiation  of  the  target  at  resonance.  There  is  a  similar  reflection  coming  from  the  bottom  of 
the  tank,  but  its  amplitude  is  considerably  smaller  due  to  the  high  absorption  of  the  redwood. 
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FIGURE  4-6.  REFLECTIONS  FROM  FRONT  WALL 

Based  on  considerations  outlined  in  the  last  few  paragraphs,  it  seems  that  a  sampling 
period  devoid  of  serious  wave  contamination  of  an  intended  signal  falls  between  the  arrival  of 
the  backscattered  returri  at  2796  )is  to  the  arrival  of  the  free  surface  reflection  at  about  4300  ^s. 
The  preferred  sampling  at  the  tail  portion  of  the  reflected  pulse  at  about  3692  \is  (see  Rgure  4-4 
(ii) )  considerably  relaxes  the  choice  of  an  allowable  transmitted  pulse  width.  As  long  as  the 
distance  from  the  sampling  location  of  the  reflected  pulse  to  its  trailing  sinusoid  is  less  than 
twice  the  target-hydrophone  separation,  no  interference  or  overlap  by  the  direct  signal  can  take 
place  irrespective  of  the  transmitted  pulse  width.  This  is  an  important  revelation  since 
conceivably  a  pulse  of  longer  pulse  width  is  needed  to  excite  large  targets  at  resonance,  and 
this  is  now  permissible.  For  most  of  the  experiments  presented  in  this  report,  however,  the 
transmitted  pulse  width  of  1  ms  seemed  adequate,  and  It  was  chosen.  This  choice  gave  an 
additional  advantage  since  the  sampling  of  the  incident  and  the  reflected  pulses  occurred  at 
different  time  periods,  allowing  each  to  be  seen  in  its  entirety  in  the  scope  screen  as  shown  in 
Figure  4-7. 

The  standard  arrangement  to  obtain  the  spectrograms  of  the  target  is  when  the  projector, 
the  hydrophone  and  the  target  are  lined  up  along  a  diameter  of  the  tank.  To  measure  the 
bistatic  scattering  cross  section  or  the  reradiation  rosetta  and  present  the  result  in  an  angular 
plot,  the  hydrophone  is  located  at  offset  positions  (see  Figure  3-6).  The  above  analysis, 
however,  can  be  applied  in  a  similar  manner.  It  can  be  shown  that  under  certain  hydrophone 
orientations,  interference  with  the  reradiation  by  the  residual  of  the  incident  pulse  is 
unavoidable.  This  will  be  addressed  later  when  the  measurement  results  are  discussed. 

The  sampling  location  in  the  reflected  pulse  depends  on  receive  time  delay  and  the 
received  pulse  width.  As  discussed  earlier,  the  first  and  last  several  cycles  are  avoided  and 
since  the  amplitude  of  the  reradiation  decreased  from  cycle  to  cycle,  it  was  necessary  to  make 
the  received  pulse  width  narrow,  yet  wide  enough  to  enhance  the  signal  to  noise  ratio.  A 
choice  of  0.1  ms  for  the  received  pulse  width  was  found  to  be  a  good  compromise  for  most  of 
the  experiments.  The  received  sampling  gate  corresponding  to  this  choice  can  be  seen  in 
lower  trace  of  Figure  4-7. 
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FIGURE  4-7.  ACTUAL  OSCILLOSCOPE  TRACE 

In  addition  to  the  received  pulse  width,  the  intensity  of  the  incident  wave  also  affects  the 
signal-to-noise  ratio,  in  general,  the  amplitude  of  the  Incident  wave  was  set,  within  a  specific 
frequency  range,  as  high  as  possible  without  incurring  saturation  and  distortion  in  the  projector 
voltage.  The  same  concept  was  applied  to  the  receiving  side  with  the  amplification  of  the 
signal  detected  by  the  hydrophone.  The  product  of  the  frequency  and  the  pulse  width  gives 
the  number  of  cycles  impinging  on  the  target  per  incident  pulse.  It  was  chosen  such  that  a 
target  always  achieved  equilibrium  by  the  end  of  the  insonifying.  A  large  target  might  need  a 
pulse  of  large  pulse  width  and  reduced  repetition  rate.  For  the  samples  used  in  the 
experiments,  a  repetition  rate  of  10  pulses  per  second  was  found  to  be  adequate. 

The  above  discussion  presents,  in  fair  detail,  a  comprehensive  description  of  the 
approach  and  procedure  in  the  choice  of  pertinent  parameters  for  an  exercise  to  extract  the 
signature  of  an  underwater  target  by  an  active  sonar  scheme  based  on  resonances.  Some  of 
the  findings,  such  as  the  relaxation  of  the  restriction  on  the  transmitted  pulse  width  and  the 
practicality  of  the  parameter  selections  based  on  the  fruitful  execution  of  the  scheme,  have 
enhanced  the  probability  of  success  with  realistic  targets  in  the  future.  Based  on  simple 
acoustical  scaling,  a  comparison  of  several  operational  parameters  of  a  typical  target  used  in 
the  present  investigation  with  two  future  candidates  is  given  in  Table  4-1 .  As  can  be  seen  from 
the  table,  it  is  possible,  without  modifying  the  present  experimental  scheme,  to  increase  the 
pulse  width  or  the  pulse  repetition  rate  by,  at  least,  a  combined  factor  of  25  for  the  36"  model. 
Under  this  new  condition,  there  will  be  some  900  sinusoidal  impingements  per  second  on  the 
target  at  2,778  Hz.  For  the  360"  model,  the  impingements  on  the  target  could  be  increased  to 
90  per  second  at  278  Hz. 
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There  is  no  restriction  on  further  increase  of  the  incident  puise  width  until  several 
resonances  of  a  target  are  secured  by  a  corresponding  decrease  in  pulse  repetition  rate.  As 
long  as  the  body  of  water  is  greater  than  X/4  in  its  least  dimension.  For  the  287  Hz  case,  this  is 
80’. 


TABLE  4-1.  PERTINENT  OPERATIONAL  PARAMETERS  IN  THE  RST  EXPERIMENTS 


Operational 

Parameters 

Present 

Model 

Future 

Model  1 

Future 

Model  2 

Size  (inches) 

1 

36 

360 

Incident  Frequency  (Hz) 

100,000 

2,778 

278 

Pulse  width  (ms) 

1 

36 

360 

Cycles  in  the  puise 

100 

100 

100 

Pulse  rate  (Hz) 

10 

10 

10 
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CHAPTER  5 

EXPERIMENTAL  RESULTS 


Many  results  were  obtained  based  on  the  wide  variety  of  targets  measured;  only  a  few 
sample  results  are  given  here.  Even  though  these  samples  obviously  are  the  best  and  the  most 
representative,  many  others  with  good  quality  exist  that  could  equally  have  been  selected. 

Most  of  the  samples  were  chosen  particularly  with  regard  to  the  verification  of  the  theory  and 
establishing  a  workable  methodology  for  the  classification  of  underwater  targets  via  active 
sonar,  which  is  the  goal  of  the  present  investigation. 

It  is  advantageous,  at  this  point,  to  took  at  some  additional  theoretical  solutions  to  be 
compared  to  the  experimental  results.  In  Chapter  2,  theoretical  plots  were  presented  with 
discussion  pertaining  to  the  subtraction  of  the  rigid  response  from  the  total  response  to  leave 
the  modal  or  resonant  responses  of  the  target.  For  simplicity,  since  there  is  only  one  k  in  the 
text  from  here  to  the  end  the  subscript  is  dropped  (k,  =k).  Rgure  5-1  shows  the  result  of 

calculations  by  Gaunaurd  and  BrilP^-'^  and  presents  the  modal  responses,  in  the  ka  domain,  of 
an  infinite  solid  cylinders  and  a  hollow  cylinder  with  aspect  ratio  of  ID/OD  =  2/3.  The  numbers 
beside  the  resonance  peaks  are  the  /  index  of  the  (n./)  mode. 

The  plots  are  for  insonification  by  a  plane  wave,  of  constant  amplitude,  broadside  to  the 
axis  of  symmetry.  It  Is  an  alternative  form  of  the  spectrogram  of  the  reradiated  return  shown  as 
the  bottom  trace  In  Figure  5-2.  The  theoretical  data  in  Figure  5-1  does  not  include  noise  or  any 
uncertainty  that  is  associated  with  the  real  world.  By  comparison  of  the  features  between  the 
theoretical.  Figure  5-1  and  the  experimental  Rgure  5-2,  such  as  the  locations  of  the 
resonances,  their  respective  amplitude  (size)  and  shape  will  determine  the  validity  of  the 
measurements  scheme  and  assumptions  made.  Note  that  the  increasing  frequency  scale  runs 
in  opposite  directions  between  the  theoretical  and  the  experimental. 

A  close  examination  of  the  theoretical  curves  In  Rgure  5-1  reveals  that  each  mode  has  its 
own  unique  shape,  amplitude,  and  location  in  the  ka  domain.  Some  resonances  are  wider  than 
others,  and  the  distance  in  the  frequency  domain  between  the  modes  is  larger  in  hollow  than  in 
solid  cylinders.  This  makes  a  hollow  tube  a  better  target  for  verification  of  the  resonance  theory. 
These  unique  characteristics  will  be  used  in  this  study  to  determine  which  mode  has  been 
experimentally  determined  and  to  compare  how  close  it  compares  to  theory.  A  review  of 
classic  vibration  theory  will  reveal  that  the  width  of  the  resonances  are  dependent  on  the  loss 
mechanisms  in  the  structure.  The  targets  used  in  this  study  have  at  least  two  mechanisms 
associated  with  the  dissipation  or  loss  of  energy,  namely  the  flow  of  energy  from  the  target 
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(radiation)  and  the  internai  energy  (heat).  The  heat  wiii  be  dissipated  to  the  surrounding  water, 
but  the  process  is  much  slower  and  is  ignored.  If  a  certain  mode  is  an  efficient  radiator,  the 
time  domain  signal  of  the  reflected  and  the  reradiated  energy  will  show  a  quick  decay  and  the 
frequency  domain  signal  will  have  a  large  or  wide  resonance  peak.  The  wide  resonance 
indicates  a  low  "Q"  or  lossy  system.  For  these  cases,  it  has  been  experimentally  determined 
that  the  input  energy  used  in  this  study  was  not  sufficient,  in  many  cases,  to  produce  a  peak 
much  above  the  noise  floor  of  the  measuring  system.  If  the  heat  is  the  dominant  term,  as  in 
inefficient  radiating  modes,  the  time  domain  signal  will  have  a  slow  decay  and  the  frequency 
domain  signal  will  have  a  very  narrow  modal  response  (high  Q  systems).  These  are  the  most 
likely  modes  to  be  observed  experimentally,  especially  the  ones  in  Figure  5-1  with  the  higher 
amplitudes.  In  these  instances,  both  the  frequency  location  and  the  other  features,  size  and 
shape,  are  expected  to  have  good  theoretical  to  experimental  correlation. 


mm 

1 

smui 

HUH! 

1 

jSjl^ 

IRR 

Hi 

HH 

1 

liSH 

IRHrI 

!■■■■■■ 

mm 

wm 

1 

■BB 

mmm 

^hh^ihhh 

mm 

mm 

1 

^HHHi^HH 

1 

IRRI 

FIGURE  5-1 .  PREDICTED  RESPONSE  IN  THE  ka  DOMAIN 


At  some  frequencies  or  ka  values,  more  than  one  mode  may  be  excited.  Not  only  would 
the  Incident  energy  in  these  cases  be  partitioned  amongst  all  the  modes  near  or  at  that  ka 
value,  but  also  the  resulting  peaks,  due  to  reradlation  from  many  modes,  may  not  resemble  any 
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one  of  the  single  theoretical  modal  responses.  This  phenomenon  occurred  ratfier  regularly  and 
might  cast  doubts  about  the  verification  of  the  RSI  predictions.  However,  there  are  remedies  in 
the  measurement  of  the  monostatic  reflections  which  can  serve  to  remove  most  of  the 
uncertainties.  As  examples  of  the  discussions  just  presented,  the  (2,1),  (2,4),  and  the  (2,7), 
modes  in  Figure  5-1  for  solid  cylinders  have  wider  (more  lossy  or  lower  Q)  than  the  (2,3)  and 
(2.5)  modes  and  are  less  likely  to  be  excited  in  an  experiment  than  the  (2,3)  and  the  (2,5) 
modes.  The  (2,2)  mode,  which  has  a  very  small  amplitude,  is  also  unlikely  to  be  observed 
because  of  the  noise  in  the  experimental  system. 


Other  calculations,  which  appear  in  the  literature,  are  much  less  popular  than  those  of  the 
modal  response  in  Figure  5-1 .  The  most  common  presentation  is  a  plot  of  the  monostatic  cross 
section  which  expresses  the  steady  state  (equilibrium)  backscattered  energy  return.  These 
calculations  predict  the  drop  in  target  strength  and  can  be  compared  to  the  middle  curve  of  the 
measured  data  plots  in  Figure  5-2.  For  solid  cylinders,  the  response  is  predicted  by  Equation 
(2-33): 


na 


£  f„(x,Tc) 


17IX»=0 


u-nw 


(2-33) 


It  is  important  to  note  that  both  the  reflected  and  the  reradiated  energy  are  affected  at  or  near  a 
resonant  frequency.  Even  though  either  can  be  used  to  identify  the  resonance,  in  this  work 
both  are  measured  and  presented  in  order  to  provide  direct  comparison  to  many  published 
predictions. 
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If  the  target  is  suspended  horizontally,  the  reflected  energy,  in  polar  form,  will  result  in  a 
pattern  that  has  broad  lobes  at  and  near  broadside  and  spikes  at  the  flat  ends.  It  is  commonly 
referred  to  as  a  "butterfly*  pattern.  At  a  resonance,  this  pattern  is  characterized  by  a  deep 
notch  in  the  backscattered  direction  at  6  =  n.  For  example,  see  Figure  Sin  the  work  by  Brill  and 
Gaunaurd'^as  well  as  Figure  5-3  in  this  report  shows  a  notch  in  the  backscattered  signal  which, 
of  course,  signifies  reduced  energy  return  at  resonance  and  a  subsequent  increase  in  the 
reradiated  energy  at  broadside.  This  leaves  little  doubt  that  the  resonances  of  cylinders  is  a 
broadside  phenomenon. 

For  a  cylinder  suspended  vertically,  the  polar  pattern,  or  directivity  pattern,  of  the 
reradiation  is  theoretically  given  by  |cosa6|  which  results  in  the  familiar  "figure  8*  pattern  for 
n  =  1  and  the  "four-leaf  clover"  pattern  for  n  =  2,  etc.  Figure  5-4  is  an  example  of  a  measured 
ros 


Many  ka  domain  or  spectral  plots  similar  in  form  to  Rgure  5-1  are  given  in  following 
figures.  Figures  5-5  pertains  to  the  solid  aluminum  cylinders  of  diameter  of  1.867*.  As  with  all 
spectral  plots,  the  range  of  the  experiment  is  given  by  the  frequency  and  corresponding  ka 
scale. 


FIGURE  5-3.  REFLECTED  AND  RERADIATED  ENERGY  VERSUS  ASPECT  ANGLE  AT 


RESONANCE 

Generally  the  experimental  results  checked  very  well  with  the  theoretical  predictions 
presented  in  Figures  5-1  for  cylindrical  targets.  Table  5-1  provides  a  tabular  comparison  of  the 
theory  and  experimental  results. 
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FIGURE  5-5.  SPECTROGRAM  FOR  A  ^.86r  BY  6"  SOLID  CYLINDER 


5-5 


css  TR  456-93 


TABLE  5-1 .  RESONANCES  OF  A  1 .876*  BY  6“  SOLID  ALUMINUM  CYLINDERS 


1  Mode  Number 

1.2 

1.3 

1,4 

MEM 

mm 

ESI 

Ea 

6,1 

ka 

r 

13.9 

16.0 

16.8 

16.2 

19.5 

ka 

M* 

14.0 

15.8 

16.9 

16.4 

19.7 

Aka 

0.1 

0.2 

0.1 

EM 

M’ 

154 

174 

187 

181 

218 

1  T-  Theoretical,  M'-  Measured  I 

As  shown  in  the  previous  figure  and  table,  the  prominent  peaks  are  at  154, 174, 187, 181 , 
and  218  kHz  correspond  to  the  (1 ,3),  (1 ,4),  (2,3),  (4,2),  and  (5,2)  mode  which  occur  at  ka=  14.0, 
15.8, 16.9, 16.4,  and  19.7  respectively.  The  variance  in  the  data  is  directly  related  to  the  width 
of  the  theoretical  resonances  as  shown  in  figure  5-1  due  to  the  low  signal  as  discussed 
previously.  The  noisiness  of  the  data  makes  determination  of  the  exact  peak  difficult. 


FIGURE  5-6.  EXPANDED  VIEW  OF  (5,2)  AND  (3,3)  MODES 

The  peak  at  154  kHz  is  due  to  the  clustering  of  the  (1 ,3)  mode  at  ka=13.9  and  the  (6.1 ) 
mode  at  ka=14.1 .  When  clustering  like  this  occurs,  the  resultant  peak  is  also  affected  by  the 
relative  amplitude  of  the  resonance  peak  in  the  cluster.  A  case  in  point  is  the  (3,3)  mode  which 
has  aka  value  of  19.6  and  is  very  near  that  for  the  (5,2)  mode.  Even  though  both  these 
resonances  are  very  narrow  and  should  have  good  experimental  results  the  combination  results 
in  a  larger  error.  The  two  resonances  are  plotted  in  Figure  5-6  with  an  expanded  scale  to  show 
this  relationship. 
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In  Figure  5-5  The  resonance  was  labeled  for  the  (5,2)  mode  since  it  is  dominant.  Figure 
5-7  is  the  spectrogram  for  a  0.978*  by  6*  solid  aluminum  cylinder  and  Table  5-2  is  the  tabular 
form. 


FIGURE  5-7.  SPECTROGRAM  FOR  A  0.978*  BY  6*  SOLID  CYLINDER 

The  prominent  resonance  frequencies  occurred  at  264, 300, 310,  and  319  kHz.  These 
frequencies  correspond  to  the  (6,1),  (1,4),  (4,2),  and  (2,3)  modes  atka=  14.0, 15.9, 16.4,  and 
16.9,  respectively.  These  values  are  within  0.1  of  theoretical.  Due  to  the  smaller  size,  the  (6,1 ) 
mode  (ka=  14.0)  and  the  (1 ,3)  mode  (ka«13.6)  can  be  determined  separately.  The  apparent 
resonance  at  288  kHz  (ka=15.2)  seemed  to  indicate  a  resonance  that  is  not  covered  by  theory. 
The  peak  at  343  kHz  (ka=18.3)  is  a  classic  example  of  a  wide  bandwidth  resonance  that  is 
difficult  to  excite.  Again,  the  flatness  and  the  noisiness  of  the  bottom  trace  in  Figure  5-7 
indicates  that  the  data  was  taken  near  the  noise  floor  of  the  measurement  equipment  and  in 
several  cases  the  reradiated  energy  was  barely  detectable  above  this  noise  floor. 

TABLE  5-2.  RESONANCES  OF  A  0.9.876"  BY  6"  SOLID  ALUMINUM  CYLINDERS 


1  Mode  Number 

1.2 

1.3 

1.4 

■9 

lii 

wa 

6,1  I 

ka 

a 

16.0 

16.2 

16.3 

14.1  1 

ka 

M* 

15.9 

16.0 

16.4 

14.0 

Aka 

0.1 

0.2 

0.1 

0.1 

Wm&lm 

M’ 

300 

319 

310 

264 

1  T=  Theoretical,  M**  Measured 

1 
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Figure  5-8  is  the  spectrogram  for  a  0.740"  by  6*  solid  aluminum  cylinder  spanning 
frequencies  from  350  to  450  kHz.  Despite  the  raised  noise  floor  on  the  data  plot  (reduced 
signal  level),  prominent  peaks  can  be  detected  at  396, 413,  and  425  kHz  which  corresponds  to 
(1,4),  (4,2),  and  (2,3)  modes  at  ka  values  of  15.8, 16.4,  and  16.8.  The  reduced  size  of  the  target 
provides  an  advantage  of  separating  the  modes  and  a  disadvantage  of  a  reduced  signal  level. 
Table  5-3  provides  a  tabular  form  of  the  data. 


FIGURE  5-8.  SPECTROGRAM  FOR  A  0.740"  BY  6"  SOLID  CYLINDER 
TABLE  5-3.  RESONANCES  OF  A  0.74"  BY  6"  SOLID  ALUMINUM  CYLINDERS 


Mode  Number 

1.2 

1.3 

1.4 

MB 

EH 

IBM 

~6.n 

ka 

n 

16.0 

16.8 

16.2 

ka 

M' 

15.8 

16.8 

16.3 

_ 

Aka 

0.2 

0.0 

0.1 

■IBH 

M' 

398 

425 

413 

HHHIIIi 

-  - - - -i 

T  =  Theoretical,  M'=  Measured 


This  concludes  the  data  for  a  target  of  a  single  geometry,  an  infinite  solid  cylinder,  and 
the  results  are  presented  in  Table  5-4  which  shows  that  indeed  the  resonances  occur  at  the 
same  ka  and  at  different  frequencies  for  the  different  diameter  targets. 
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TABLE  5-4.  RESONANCES  OF  SOLID  CYLINDERS  WITH  THREE  DIFFERENT  DIAMETERS 


Figure  5-9  shows  the  comparison  of  the  theoretical  and  the  measured  shape  of  the  (1.2) 
and  the  (2,3)  mode  of  a  0.989"  diameter  solid  aluminum  cylinder.  The  agreement  is  quite  good. 


FIGURE  5-9.  DETAILED  SPECTROGRAM  OF  A  (1 .2)  AND  A  (2.3)  RESONANCE 


Figure  5-10  comprises  three  spectral  plots  for  a  hollow  aluminum  cylinder  with  an  ID/OD 
ratio  of  2/3  and  a  diameter  of  0.75".  The  three  plots  cover  the  frequency  ranges  from  50  kHz  to 
450  kHz  or  ka  range  from  2  to  18.  The  convergence  and  the  reduction  in  amplitude  of  the 
direct,  reflected  and  reradiated  traces  sXka=2  is  an  indication  that  we  are  reaching  the  limit  of 
the  test  set-up  and  higher  power  levels  are  required. 
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FIGURE  5-10.  SPECTROGRAM  OF  A  0.75*  DIAMETER  HOLLOW  CYLINDER 
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Table  5-5  is  a  tabular  form  of  the  data  from  the  previous  figure.  Figure  5-1 1  provides  the 
spectrogram  for  a  1 .5*  diameter  hollow  cylinder  with  the  same  aspect  ratio  (2/3)  over  a  range  of 
frequencies  from  50  to  400  kHz  and  Table  5-6  provides  the  tabular  form  of  the  data.  The 
spectral  plots  for  the  hollow  tubes  show  more  widely  spread  peaks  and  dips  as  predicted  by 
the  theory.  Table  5-7  provides  a  comparison  between  the  two  cylinders  previously  presented. 
Again  the  ka  values  agree  in  each  mode  and  the  frequency  is  different.  Measurements  for 
these  two  cylinders  revealed  only  three  common  modes. 

TABLE  5-5.  RESONANCES  OF  A  0.75"  DIAMETER  HOLLOW  ALUMINUM  CYLINDER 


1.4  I  2.2  I  2.3  3.2  3.3  1.3  5.2  6.1 


T  =  Theoretical.  M‘=  Measured _ 


TABLE  5-6.  RESONANCES  OF  A  1 .5"  DIAMETER  HOLLOW  ALUMINUM  CYLINDER 


2.3 

3.3 

4.2 

4,3 

5.2 

5,3 

23.0 

25.4 

16.4 

28.6 

20.0 

31.4 

23.1 

25.6 

16.6 

28.6 

20.0 

31.4 

0.1 

0.2 

0.2 

0.0 

0.0 

0.0 

286 

318 

208 

355 

258 

390 

T  =  Theoretical.  M’=  Measured  _ 


TABLE  5-7.  RESONANCES  OF  HOLLOW  CYLINDERS  WITH  TWO  DIFFERENT  DIAMETERS 


6.1 


Mode  Number 


Frequency  1^ 

(kHz)  2* 


ka 


234  322  412 


115  286  318  208  |  355  |  248  |  390 


16.6 


25.6  I  16.6  28.6  I  20.0  |  31.4 


1*=>  Diameter=0.75".  2‘=>  Diameter=1.5" 
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FIGURE  5-1 1 .  SPECTROGRAM  OF  A  VERTICAL  1.5"  DIAMETER  HOLLOW  CYLINDER 
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The  plots  up  to  this  point  pertain  to  targets  suspended  with  their  axis  in  the  vertical 
position.  In  Figures  5-12  the  same  1 .5“  hollow  cylinder  with  aspect  ratio  of  2/3  was  suspended 
with  its  axis  in  a  horizontal  position.  Since  the  case  of  normal  incidence  prevailed  they  illustrate 
the  effect,  or  the  lack  of  it.  of  the  positioning  of  the  target. 


FIGURE  5-12.  SPECTROGRAM  OF  A  HORIZONTAL  1.5"  DIAMETER  HOLLOW  CYLINDER 


Next  supporting  data  will  be  supplied  to  prove  that  the  circumferential  directivity  pattern 
of  the  reradiation  is  given  by  |  cosn0|  which  results  in  the  familiar  "figure  8"  pattern  for  n  =  1  and 
the  "four-leaf  clover"  pattern  for  n  =  2,  etc.  Figure  5-13  shows  two  "Rosetta"  patterns  for  the 
same  cylinder  at  two  different  frequencies.  The  theoretical  curve  is  provided  on  portions  of  the 
curve  In  dotted  lines.  The  oscillations  around  0  =  0®  is  due  to  interference  between  the  residual 
of  the  direct  energy  and  the  reradiated  energy. 

Figure  5-14  is  a  comparison  of  the  "Rosettas"  for  two  different  diameters  (1.5"  and  6.0") 
hollow  aluminum  cylinders  with  the  same  aspect  ratio  of  2/3. 
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(4.2)  Mode  (4,3)  Mode 

FIGURE  5-13.  ROSETTAS  OF  THE  (4.2)  AND  (4,3)  MODE  FOR  A  1.5"  HOLLOW  CYLINDER 


FIGURE  5-14.  ROSETTAS  OF  THE  (5.2)  MODE  FOR  TWO  DIFFERENT  DIAMETERS 


Figure  5-15  is  the  "Rosettes"  for  the  (2,3)  and  (3,3)  mode  for  1 .5"  hollow  aluminum  cylinder 
with  aspect  ratio  of  2/3  showing  that  the  first  index,  n,  is  for  the  circumferential  wave  which 
affects  the  rosetta  while  the  second  index,  C,  is  for  the  internal  waves  which  is  not  propagated 
into  the  water. 


css  TR  466-93 


FIGURE  5-15.  ROS^AS  OF  THE  (2.3)  AND  (3.3)  MODE  FOR  A  HOLLOW  CYLINDER 


As  stated  previously  some  times  there  are  more  than  one  mode  present  at  a  given  ka 
location  and  the  resultant  shape  and  directivity  pattern  will  be  the  sum  of  all  modes  excited. 
Figure  6-16  is  a  good  example  of  that  process.  The  measured  data  is  given  in  the  solid  lines 
and  the  theoretical  (in  the  upper  left  hand  quadrant)  is  given  as  a  dashed  line  for  a  1 .5" 
diameter  hollow  aluminum  cylinder  with  aspect  ratio  of  2/3  measured  at  161kHz.  The 
agreement  is  good. 

The  final  experiment  associated  with  cylinders  was  to  measui’e  the  effect  of  the  target 
being  close  to  but  not  at  broadside.  Figure  5-17  presents  the  results  of  a  very  narrow  portion  of 
the  spectrogram  of  a  1.5"  hollow  cylinder  with  the  (2,3)  mode  excited.  The  target  was  tested  at 
0®.2‘’  ,4®,  and  10®  from  broadside.  The  trend  of  the  data  was  for  the  resonance  to  move  to  a 
lower  frequency  and  for  its  amplitude  to  decrease  until  at  10®  the  resonance  is  barely 
distinguishable. 
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FIGURE  5-17.  DEPENDENCE  OF  THE  (2.3)  MODE  ON  ANGLE  OF  INCIDENCE 
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Considerations  and  measurements  to  this  point  have  been  associated  with  cylinders 
which  are  strictly  a  broadside  phenomen4Di.This  is  not  the  case  in  the  Figures  5-18  through 
5-22,  which  provides  spectral  results  from  nrK)re  realistic  targets:  a  1 .4“  by  T  solid  spheroid  and 
a  3.5“  by  7"  hollow  spheroid  with  their  major  axes  in  the  horizontal  orientation.  Figure  5-18 
presents  the  backscattered  curve  and  the  spectrogram  for  the  1 .4*  by  7“  solid  aluminum 
spheroid  at  normal  incidence  to  its  major  axis  (broadside),  while  Figure  5-19  shows  the  same 
when  wave  incidence  is  at  the  nose  of  the  spheroid. 


FIGURE  5-18.  SPECTROGRAM  OF  A  1.4"  BY  T  SOLID  SPHEROID  AT  BROADSIDE 
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The  next  three  figures  pertain  to  spectral  results  from  a  3.6*  by  7“  hollow  brass  ellipsoid. 
Broadside  impingement  was  in  the  normal  direction  for  the  results  in  Figure  5-20  and  2.5®  and 
10®  off  the  normal  for  those  in  Figures  5-21  and  5-22.  Significant  matches  of  peaks  and  dips 
could  be  seen  in  all  the  figures.  The  results  further  indicate  that  resonances  occur  irrespective 
of  the  target  and  its  composition.  They  also  present  a  challenge  to  the  classification  process 
because  the  unique  signature  may  well  be  orientation  dependent. 


FIGURE  5-20.  SPECTROGRAM  OR  A  3.5"  BY  7*  SPHEROID  AT  BROADSIDE 


FIGURE  5-21  .SPECTROGRAM  OF  A  3.5"  BY  7"  SPHEROID  AT  2.5®  OFF  BROADSIDE 
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FIGURE  5-22.  SPECTROGRAM  OF  A  3.5"  BY  7"  HOLLOW  SPHEROID  AT  10"  OFF 


BROADSIDE 


Figure  5-23  is  the  polar  patterns  for  the  target  strength  and  the  re-radiated  energy  from 
the  hollow  brass  spheroid  at  resonance  at  lOl.SkHz  at  a  resonance  shown  in  Rgure  5-22. 


FIGURE  5-23.  REFLECTED  AND  RERADIATED  PATTERN  FOR  A  3.5"  BY  7"  HOLLOW 
SPHEROID  AT  BROADSIDE 
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The  resonance  phenomenon  can  also  be  observed  from  the  reflected  time  domain  waves. 
As  mentioned  at  several  occasions  previously,  the  envelope  of  the  reflected  waves  virtually 
remains  rectangular  in  shape  at  nonresonant  states,  but  will  gradually  deviate  from  it  at 
resonance  due  to  the  decrease  in  amplitude  In  the  trailing  cycles.  The  decrease  signifies 
energy  penetration  into  the  elastic  target.  Moreover,  at  the  end  of  the  reflected  pulse  with  no 
significant  input  forthcoming,  the  sudden  burst  of  signal  picked  up  by  the  hydrophone  is  due  to 
the  reradiation  from  the  ringing  target.  Shown  in  Rgure  5-24  are  two  sets  of  pictures  to  illustrate 
this  phenomenon.  The  three  pictures  at  the  left  pertain  to  the  0.74*  solid  aluminum  rod  at 
frequencies  of  152, 153  and  154  kHz.  The  resonance  frequency  of  153  kHz  corresponds  to  the 
(1,2)  mode  at  ka  =  6,2.  The  three  pictures  at  the  right  pertain  to  the  1 .5"  hollow  aluminum  tube  of 
ID/OD  =  2/3  at  frequencies  of  285, 286  and  287  kHz.  The  resonance  frequency  of  286  kHz 
corresponds  to  the  (2,3)  mode  at  ka  =  22.7.  Both  sets  vividly  illustrate  the  transitions  towards 
and  away  from  the  resonance  with  the  accompanying  decrease  in  amplitude  at  the  trailing 
portion  of  the  pulse  and  the  reradiated. 


FIGURE  5-24.  REFLECTED  AND  RERADIATED  TIME  DOMAIN  RESPONSE  AT 
RESONANCE  (fj  AND  ±  1kHz. 
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Resonances  occur  at  a  set  of  ka  values  derived  from  the  solution  of  the  characteristic 
equation,  Equation  (2-53),  where  k  is  the  wave  number  and  a  is  the  radius  of  the  target.  Since 
the  wave  number  is  proportional  to  the  incident  frequency,  it  is  possible  to  achieve  a  specific 
resonance  at  various  incident  frequencies  by  using  targets  of  appropriate  sizes.  This  was 
indeed  done  in  the  experiments.  Shown  in  Rgure  5-25  are  backscattered  and  reradiated 
pulses  for  the  (4,2)  resonance  with  ka=  16.2.  The  pictures  marked  (a),  (b)  and  (c)  pertain  to 
hollow  aluminum  tubes  of  ID/OD  =  2/3;  of  outside  diameters  of  6*,  3",  and  1 .5"  respectively.  The 
corresponding  incident  frequencies  are  52.0, 102.3  and  206.9  kHz.  The  pictures  do  not  look  as 
dramatic  as  those  shown  in  Figure  5-25  but  nevertheless  exhibit  the  general  features  pertinent 
to  a  resonance  state.  The  deficiencies  often  resulted  from  tuning  problems  and  dynamic 
limitations  in  the  projector  and  the  hydrophone  at  the  resonance  frequencies  which  led  to  low 
signal  to  noise  ratio. 


FIGURE  5-25.  REFLECTED  AND  RERADIATED  TIME  DOMAIN  RESPONSE  OF  SIMULAR 


TUBES 

The  most  interesting  feature  in  the  pictures  is  the  reradiated  pulse  which  immediately 
follows  the  reflected  pulse.  An  enlarged  picture  of  this  portion  pertaining  to  Figure  5-25(c)  is 
given  in  Figure  5-25(d).  The  reradiated  displays  the  rate  of  energy  transmission  to  the 
surrounding  inviscid  fluid  which  is  governed  by  the  dissipation  mechanism  of  the  target. 
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If  the  history  of  the  energy  dissipation,  given  by  the  envelope  of  the  reradiated  decay  is 
available,  a  simple  mathematical  derivation  will  show  that  the  damping  constant,  can  be 
given  by  any  two  amplitudes  xo  and  x„,of  the  envelope,  namely: 


1 

2rcfAT 


ln(Xo/xJ 


(5-1) 


where  f  is  the  resonance  frequency  and  AT  is  the  time  lapse  for  the  amplitude  decay  from  xo  to 
X..  By  using  Equation  (5-1)  repeatedly  over  many  of  these  envelopes  at  various  resonances,  it 
is  possible  to  obtain  damping  information  over  a  wide  range  of  frequencies. 


Finally  an  experiment  can  be  conducted  to  see  if  this  resonance  could  be  exploited  to 
determine  the  damping  characteristics  of  a  system.  Two  separate  schemes  are  envisioned. 
Rrst  after  partitioning  the  energy  radiated  from  the  system  and  that  dissipated  in  the  system, 
the  damping  factor  of  the  material  can  be  determined.  The  second  scheme  involves  selecting 
a  target  with  a  resonance  near  the  frequency  of  interest,  measuring  that  resonance  for  a  bare 
metal  target,  applying  a  coating  to  either  the  inside  or  outside  dependent  on  the  application, 
measuring  the  system  with  the  coating,  and  finally  comparing  the  results.  This  comparison  can 
be  done  in  the  time  domain  by  measuring  the  exponential  decay  of  the  reradiated  energy  of  in 
the  frequency  domain  by  looking  at  the  change  in  width  of  the  resonance  at  the  3-dB  down 
point. 


The  portion  of  the  time  domain  signal  pertaining  to  the  free  ringing  can  be  used  to 
determine  the  damping  constant  of  the  target  material.  A  lightly  damped  exponential  decay 
curve  is  presented  in  Figure  5-26.  Several  points  on  the  curve  are  plotted  in  a  semilog  plot  in 
Figure  5-27.  The  result  of  the  plot  yielded  a  damping  ratio  of  0.072  by  Equation  (5-1 ).  This  is 
probably  a  reasonable  good  value  for  aluminum  but  the  radiation  impedance  of  the  process 
needs  to  be  established. 
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FIGURE  5-26.  LIGHTLY  DAMPED  RERADIATED  PULSE 


FIGURE  5-27.  SIMILOG  PLOT  OF  LIGHTLY  DAMPED  RERADIATED  PULSE 
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CHAPTER  6 

CONCLUSIONS  AND  RECOMMENDATIONS 


In  conclusion,  the  theoretical  responses  of  simple  shapes  were  compared  to  the 
measured  values  and  the  agreement  was  found  to  be  excellent.  More  complex  structures  like 
spheroids  were  measured  and  were  also  found  to  have  detectable  resonances.  A  cursory 
review  of  the  frequency  dependence  of  the  resonances  at  angles  other  than  broadside  and  a 
scheme  to  use  the  resonances  to  determine  the  damping  characteristics  of  the  targets  was 
completed.  This  body  of  experiments  has  proven  beyond  a  shadow  of  doubt  that  the 
resonances  of  simple  structures  can  be  mathematically  determined  and  that  they  can  be 
measured.  Resonances  were  measured  for  solid  and  hollow  cylinders  as  well  as  for  spheroids. 
The  latter  data  were  presented  without  analysis.  In  summary,  the  following  extension  of  the 
present  experiments  is  recommended: 

The  resonances  of  more  complex  structures  need  to  be  measured.  This  would  involve 
building  several  small  targets  which  resemble  a  target  of  interest  with  varied  size.  The 
experimentally  determined  set  of  resonances  would  become  the  "class*  signature  of  that  target. 

To  investigate  the  application  of  this  scheme  for  identification  of  mines,  a  mine  model 
should  be  built  and  tested  in  the  free  field  and  on  the  bottom  to  determine  how  the  bottom 
interaction  affects  the  resonance.  It  is  suspected  that  an  increase  in  the  damping  constant  will 
be  noted. 

Perform  additional  experiments  and  analysis  on  the  damping  characteristic  of  metals. 

This  scheme  may  prove  valuable  as  a  means  of  measuring  the  damping  constant  of  a  fluid 
loaded  coated  cylinder. 

Finally  an  investigation  of  using  an  insonification  pulse  that  is  broadband  or  band  limited 
in  nature.  The  band  of  the  signal  should  cover  the  class  signature  of  the  target  previously 
discussed.  This  scheme  should  also  use  the  Wigner-Ville  distribution  scheme  to  analyze  the 
sonar  returns.  The  Wigner-Ville  distribution  function  is  very  easy  to  implement  and  it  determines 
the  time  dependence  of  the  frequency  domain  signal.  If  intensity  modulation  is  used  on  a  water 
fall  display,  bright  spots  will  occur  at  the  frequency  where  reflection  or  reradiation  are  present, 
and  dark  spots  other  places.  At  frequencies  associated  with  the  resonances  of  the  target  a 
bright,  dark  then  bright  sequence  will  occur  in  a  line  across  the  display.  If  this  matches  the 
class  signature  then  a  classification  can  be  made.  The  reverberation  will  show  up  as  a  solid 
bright  line  for  all  frequencies. 
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The  Navy  has  invested  heavily  in  the  task  of  underwater  target  classification  and 
identification,  and  it  seemes  that  this  trend  will  be  continued  for  years  to  come.  The 
methodology  herein  developed  is  a  very  small  part  of  the  overall  endeavor  and,  in  the  opinion  of 
the  authors,  has  shown  promise  of  fruitful  applications  in  real  situations.  It  deserves  close 
scrutiny  for  its  practicability  and  potential  for  future  development. 

Finally  Table  6-1  presents  all  the  resonances  that  were  measured  during  this  research 
and  reported  in  this  report.  Many  other  resonance  modes  have  been  determined.  "H“  signifies 
a  hollow  cylinder  and  *S"  a  solid  cylinder. 

TABLE  6-1.  EXPERIMENTALLY  DETERMINED  RESONANCES  MODES 


1  Mode  #  1 

1.2 

1.3 

1.4 

1^ 

3,2 

KSi 

KEI 

5,2 

5.3 

6.1 

H,S 

S 

S 

H,S 

H,S 

H 

H,S 

H,S 

H 

H,S 

H 

S 
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